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1.0  Overview 


The  Space  Data  Analysis  Laboratory  (SDAL)  of  Boston 
College  was  contracted  by  the  Air  Force  Geophysics  Laboratory 
(AFGL)  Aerospace  Engineering  Division  (Data  Systems  Branch) 
to  perform  multiple  tasks  related  to  spacecraft  missions. 
These  tasks  included  the  developing  of  data  bases  for  a 
multitude  of  scientific  instruments;  developing  and 
implementing  analyses  using  the  data  bases  as  the  prime 
input;  and  providing  significant  input  to  Data  Management 
Plans  for  several  missions. 

The  contract  number  under  which  these  efforts  were 
performed  was  F19628-83-C0100 . 

Data  systems  were  designed  for  several  satellite 
projects.  Satellites  involved  in  these  efforts  include  HILAT, 
Polar  Bear,  and  S85-1.  Systems  were  also  designed,  developed 
and  implemented  for  flights  of  the  Space  Transportation 
System.  In  addition,  analytic  systems  were  setup  for  several 
rocket  projects. 

System  design  used  modular  programming  techniques  to 
provide  the  necessary  flexibility  in  satisfying  requirements 
in  a  cost  effective  manner.  Subroutine  libraries  were 
developed  for  use  in  the  overall  systems.  These  libraries 
allow  for  the  performance  of  common  tasks  and  assure  that 
these  tasks  are  executed  in  identical  manners. 


2.0  HILAT  Analyses 


The  HILAT  spacecraft  designated  P83-1  was  launched  into 
polar  orbit  on  a  Scout  rocket  on  27  June  1983  from  Vandenburg 
AFB .  The  HILAT  payload  was  designed  to  provide  data  for  the 
study  of  plasma  density  irregularity  dynamics  in  the  high 
latitude  ionosphere.  The  spacecraft  orbit  is  near  circular 
at  800km  with  an  orbital  period  of  approximately  101  minutes. 
In  addition,  the  vehicle  is  near  three  axis  stabilized 
(drifts  in  pitch,  yaw  and  roll  can  be  seen  on  each 
acquisition) . 

HILAT  operates  in  real  time  mode  only  with  data 
transmitted  to  any  o t  four  ground  stations.  Three  of  these 
stations  are  at  fixed  high  latitude  sites  (Tromso,  Norway; 
Sondrestrom,  Sweden;  and  Fort  Churchill,  Canada) .  The 
remaining  station,  refered  to  as  the  Rover,  is  moveable  from 
location  to  location  although  it's  most  frequent  site  is 
Seattle,  Washington. 

The  HILAT  payload  consists  of  a  coherent  beacon  which 
provides  scintillation  and  total  electron  content  data;  a 
Thermal  Plasma  Monitor  which  consists  of  an  Electron  Sensor, 
Retarding  Potential  Analyzer  (RPA) ,  and  Ion  Driftmeter  ( IDM) ; 
an  Electron  Spectrometer  for  the  measurement  of  electron  flux 
in  the  energy  range  between  20  eV  and  20  keV;  the  Auroral 
Ionospheric  Mapper  (AIM)  which  consists  of  a  scanning  imager 
(which  failed  on  23  July  1983)  and  two  photometers;  and  a 
triaxial  fluxgate  magnetometer.  Attitude  data  is  computed 
using  the  triaxial  magnetometer  along  with  a  sun  sensor. 

For  each  station  pass,  the  downlinked  data  is  digitized 
by  the  station  computer  onto  a  disk  file  of  the  raw  telemetry 
data.  Between  station  passes,  the  raw  data  is  semi-processed 
and  put  onto  a  disk  file  in  "Summary  Tape"  form.  As  time 
permits,  the  raw  data  files  are  off-loaded  onto  "Raw  Tapes" 
and  the  "Summary  files"  off-loaded  onto  tapes  refered  to  as 
Summary  Tapes. 


Data  on  the  raw  tapes  consists  of  station  coordinates, 
Keplerian  elements,  spacecraft  position,  time  tagged  raw 
telemetry  data  in  packed  form,  and  raw  scintillation  data. 

Raw  tapes  from  Tromso  are  received  at  AFGL  where  they  are 
logged  and  archived. 

The  Summary  Tape  structure  has  evolved  into  a  form  of 
interlaced  records  which  contain  time,  geometry  and  model 
magnetic  field  data;  raw  telemetry  data;  scintillation 
summary  data;  and  science  summary  data.  Each  file  contains 
data  from  one  station  pass.  Summary  tapes  from  all  stations 
are  received  at  AFGL  where  they  are  quality  checked,  copied 
for  dissemination  to  other  user  agencies  and  used  as  the 
basic  input  to  the  Scientific  Satellite  Data  Analysis  System 
(SSDAS)  efforts.  SDAL  personnel  were  deeply  involved  in  the 
design  of  and  parameter  selection  for  the  Summary  tapes. 

Under  this  contract  SDAL  personnel  were  responsible  for 
the  design,  development  and  implementaion  of  the  SSDAS  for 
H1LAT .  This  system,  as  with  SSDAS  design  for  other  vehicles, 
was  built  in  modular  form  to  minimize  effort  duplication 
where  commonality  of  requirements  exist;  to  insure  that 
common  efforts  are  performed  in  exactly  the  same  way  in  all 
routines;  and  to  provide  the  flexibility  necessary  as 
spacecraft  anomalies  occur  or  requirements  change.  A  library 
of  HILAT  routines  has  been  developed  which  is  used  by  all 
processing  routines  as  necessary. 

Upon  receipt  of  a  shipment  of  Summary  tapes,  the  quality 
check  (Q/C)  and  tape  copy  effort  must  be  immediately 
initiated  since  the  tape  copies  are  to  be  shipped  to  the 
appropriate  outside  agencies  for  performance  of  their  own 
analyses.  The  Q/C-Copy  routine  was  designed  to  also  produce 
a  Summary  Tape  log  which  includes  tape  number  and  start  and 
stop  times  of  all  passes  on  each  tape.  One  such  tape  log  is 
created  for  data  from  each  station. 


Since  the  tapes  are  written  on  a  Hewlett  Packard 
computer  with  various  combinations  of  16  and  32  bit  fixed  and 
floating  point  words,  the  initial  SSDAS  effort  consists  of 
reformatting  the  summary  tapes  for  use  on  the  central  site 
computer.  Upon  completion  of  the  reformatting  effort,  the 
original  Summary  Tapes  are  archived. 

The  AFGL  payloads  on  HILAT  consist  of  the  AIM,  Electron 
Spectrometer  and  Thermal  Plasma  Monitor.  Data  bases  for  these 
instruments  along  with  the  magnetometer  and  scintillation 
data  are  created  in  time  history  form.  These  data  bases  are 
then  used  as  input  to  modular  routines  which  produce 
microfiche  displays  of  computed  parameters.  For  the  AIM, 
Electron  Spectrometer,  Thermal  Plasma  Experiment  and 
Magnetometer  two  types  of  displays  are  produced;  summary  and 
detailed  data.  The  summary  displays  consist  of  averaged  data 
for  entire  passes.  Thus,  these  plots  provide  a  survey  of  the 
instrument  data  on  a  pass  by  pass  basis.  The  detailed 
displays  contain  all  point  data  plotted  on  a  sequence  of  time 
overlapped  frames  for  each  pass.  These  displays  are  useful 
for  event  studies.  The  displays  for  each  of  the  instruments 
are  on  common  time  scales  to  allow  for  correlation. 

For  the  Thermal  Plasma  Experiment,  individual  modules 
were  developed  for  each  of  the  three  independent  experiments. 
For  these,  data  is  readout  of  telemetry  in  64  second  blocks 
which  can  be  further  subdivided  into  8  second  segments.  The 
initial  effort  consists  of  identifying  the  beginning  of  an  8 
second  block  and  determining  the  subsection  number  within  the 
64  second  cycle.  Once  found,  the  data  must  be  accumulated  in 
64  second  groups.  These  64  second  sections  contain  all 
readouts  for  the  Electron  Sensor,  RPA  and  IDM. 

The  Electron  Sensor  is  used  to  obtain  electron  density 
and  frequency  data.  For  the  swept  portion  of  the  cycle,  the 
beginning  of  each  sweep  must  be  determined.  Once  found  the 
telemetry  values  are  converted  to  sweep  voltage.  The  sensor 
output,  which  is  readout  4  times  more  frequently  than  the 
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sweep  voltage,  is  converted  to  telemetry  voltage  and  then  to 
current.  The  sensor  potential  must  be  calculated  from  the 
subcommutated  data.  A  modular  routine  developed  for  electron 
density  determination  then  receives  as  input  the  array  of 
ordered  pairs  of  sweep  voltage  and  current  along  with  sensor 
potential  and  returns  electron  density,  electron  temperature, 
and  vehicle  potential  with  respect  to  the  plasma.  Electron 
frequency  data  is  also  obtained  at  4  fixed  frequency  values. 
The  frequency  determinations  must  first  be  made.  The  root- 
mean-square  (RMS)  values  for  each  of  the  frequencies  are  then 
computed. 

The  RPA  data  consists  of  2  and  4  second  sweeps.  The 
beginning  of  each  2  and  4  second  sweep  must  first  be 
determined.  Once  determined,  the  sweep  voltage  values  are 
converted  to  actual  sweep  voltage  and  the  sensor  output 
converted  to  current.  A  modular  routine  then  receives  as 
input  the  ordered  pairs  of  voltage  and  current  along  with 
total  spacecraft  velocity  and  sensor  potential.  The  routine 
then  computes  and  returns  the  ion  temperature  and  drift 
velocity  of  ions  normal  to  the  RPA  aperture. 

The  IDM  cycle  is  8  seconds.  The  horizontal  and  vertical 
telemetry  values  must  first  be  identified  and  converted  to 
horizontal  and  vertical  drift  tangents.  A  combination  of  8,  9 
and  10  bit  words  must  be  deconvolved  from  exponent  and 
mantissa  telemetry  in  order  to  obtain  offset  and  re-zero 
values  which  are  used  in  the  calculation  of  tangent  values. 
From  the  tangent  data  another  modular  routine  is  used  which 
applies  formulae  to  compute  ion  density  and  drift  velocities. 

A  data  base  of  TPE  parameters  along  with  selected 
ephemeris  and  magnetic  parameters  is  created. 

The  Electron  Spectrometer,  refered  to  as  the  J-Sensor, 
has  three  modes  of  operation.  Three  sensors  are  mounted  at 
zenith,  40  degrees  from  zenith,  and  nadir.  Mode  1  consists 
of  zenith  data  with  8  points  per  spectra;  mode  2  also 
consists  of  zenith  data  with  16  points  per  spectra;  in  mode  3 

5 


msmimsmmsisa 


•W‘  UUWlJWlPIW^'Vllli'IV'W  '•■II*  ■  ■'JMIF'IJI  "P  '<■  •'■■■« 


data  from  all  three  detectors  is  acquired.  Bit  flags  identify 
'even'  and  'odd'  blocks.  Initially,  even  blocks  must  be 
identified  in  order  to  determine  the  instrument  mode.  Once 
identified,  data  must  be  accumulated  into  1  second  blocks 
since  this  is  the  instrument  cycle.  The  J-Sensor  data  base 
consists  of  structured  records  which  identify  instrument 
mode,  the  deconvolved  instrument  counts  arranged  from  lowest 
to  highest  energy,  decompressed  magnetometer  counts,  and  the 
required  ephemeris  parameters.  This  data  base  has  been 
stored  in  counts  because  of  the  fact  that  for  particle 
detectors,  calibrations  are  likely  to  change  as  a  function  of 
time  from  launch  and  had  calibrations  been  directly  applied 
yet  changed,  data  base  recreation  would  have  been  necessary. 

A  modular  routine  exists  which  has  the  calibration  parameters 
stored  as  a  function  of  time  from  launch.  This  routine  is 
applied  whenever  differential  flux  or  distribution  function 
computations  are  required.  For  instance,  in  the  summary  plot 
display  routine  total  flux,  total  energy  flux  and  total 
energy  are  computed  by  integrating  the  spectra.  This  display 
routine  utilizes  the  calibration  subroutine.  The  deconvolved 
magnetometer  counts  were  included  in  the  J-Sensor  data  base 
because  magnetic  pitch  angle  is  often  required  in  the 
interpretation  of  differential  flux  data.  Again,  the 
magnetometer  data  was  retained  in  counts  to  allow  for 
modified  calibrations.  A  simple  change  to  the  calibration 
module  has  taken  care  of  upgrades  to  the  magnetometer  and  J- 
sensor  calibrations. 

For  the  AIM  sensors,  individual  modules  were  developed 
to  handle  the  various  sensor  modes  (i.e.  imaging, 
spectrometer  and  photometer) .  This  proved  to  be  a  most  useful 
design  concept  since  the  imaging  mode  failed  approximately 
one  month  after  launch.  The  AIM  data  base  has  been 
succesfully  created  for  all  operational  passes.  Upon 
notification  that  the  imager  had  failed,  the  "CALLS"  to  the 
imager  subroutine  were  simply  deleted  from  the  AIM  routine  as 
opposed  to  rewriting  the  routine  to  accomodate  the  failure. 
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The  photometer  portion  of  the  routine  remains  intact  as 
originally  designed  and  the  photometer  data  base  continues  to 
be  created.  This  data  base  consists  of  the  photometer  counts 
along  with  required  ephemeris  and  attitude  parameters.  A 
calibration  subroutine  is  used  whenever  calibrated  photometer 
data  in  Rayleighs  is  required. 

The  scintillation  data  base  continues  to  be  created  from 
parameters  contained  on  the  time/geometry/magnetic  field 
records  along  with  the  scintillation  records  from  the  summary 
tape.  A  by-product  of  this  data  base  is  the  scintillation 
statistics  data  base.  This  data  base  is  created  by  inputting 
the  overall  scintillation  data  base  and  extracting  segments 
based  on  quantized  time  periods  and  averaging  the  parameters 
on  ephemeris  constraints.  Analysis  techniques  and  computer 
routines  have  been  developed  to  input  the  statistics  data 
base  and  perform  the  necessary  statistical  analyses. 

The  appendix  contains  a  sample  of  the  HILAT  tape  log, 
the  Summary  Tape  format  and  the  data  base  formats  for  the 
various  experiments. 
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3.0  Polar  Bear  Analyses 


The  Polar  Bear  satellite  was  launched  on  a  Scout  rocket 
i  in  the  fourth  quarter  of  1986.  The  spacecraft  is  very  similar 

i 

1  in  design  to  HILAT.  It  was  launched  from  Vandenburg  AFB  with 

1  an  inclination  of  nearly  90  degrees  and  an  orbital  period  of 

I  approximately  105  minutes.  Mission  duration  is  estimated  for 

3  years. 

The  requirements  for  Polar  Bear  were  similar  to  those  of 
HILAT:  Raw  and  Summary  tapes  were  to  be  archived;  data  was  to 
be  quality  checked  and  prepared  for  shipment  to  outside 
agencies;  and  a  scintillation  data  base  was  to  be  created.  In 
addition,  however,  a  data  base  was  to  be  created  for  the 
Auroral  Ionospheric  Remote  Sensor  (AIRS) .  Since  the  common 

I  requirements  between  HILAT  and  Polar  Bear  were  discussed  in 

i 

the  preceding  chapter,  discussions  here  will  concern  Polar 
.  Bear  in  general  and  AIRS  in  particular. 

j  The  telemetry  system  aboard  the  vehicle  is  similar  to 

HILAT  and  is  refered  to  as  the  Science  Data  Formatter  (SDF) . 
The  SDF  outputs  science  data  at  a  rate  of  4.098kbps.  The 
spacecraft  is  3  axis  stabilized  to  +/-  10  degrees  in  the 
j  pitch,  yaw  and  roll  axes. 

The  payload  for  Polar  Bear  includes  the  AIRS,  the  beacon 
and  the  triaxial  fluxgate  magnetometer.  The  AIRS  is  a  more 
sophisticated  version  of  the  Auroral  Ionospheric  Mapper  (AIM) 
detector  flown  on  HILAT.  It  obtains  simultaneous 
mult ispectral  UV  and  visual  images  for  both  daylight  and 
night  operations.  The  beacon,  similar  in  design  to  the  one 
flown  on  HILAT,  is  used  to  measure  scintillation  effects  at 
radio  frequencies  used  for  radar  and  communications.  The 
magnetometer  is  also  of  similar  design  to  the  system  used  on 
HILAT.  It  measures  field  aligned  currents  as  well  as 
providing  a  means  of  attitude  determination  along  with  the 
digital  sun  sensor.  These  three  instruments  will  be  operated 
concurrently  in  real  time  mode  over  the  remote  stations 


located  in  the  northern  polar  region.  The  AIRS,  however,  can 
be  operated  independently  over  mid-latitude  regions. 

There  are,  thus,  great  similarities  between  Polar  Bear 
and  HILAT.  The  Polar  Bear  Scientific  Satellite  Data  Analysis 
System  (SSDAS)  consists  of  a  modified  version  of  the  HILAT 
system.  Many  of  the  modular  routines  developed  for  HILAT 
have  direct  application  for  Polar  Bear. 

The  input  data  for  SSDAS  applications  is  the  Summary 
Tape  which  is  in  a  similar  format  to  that  used  for  HILAT. 

The  AIRS  data  base  is  being  archived  on  digital  tape.  A 
log  of  digital  tape  number  and  the  passes  contained  on  the 
tape  is  maintained.  Separate  logs  will  exist  for  each  of  the 
remote  stations  receiving  the  AIRS  data. 

The  AIRS  sensor  consists  of  four  integrated  detectors 
each  of  which  obtains  326  pixel  values  per  line  scan.  The 
time  duration  of  a  line  scan  is  2.9987  seconds  which  is 
readout  over  6  SDF  telemetry  frames.  Also  contained  in  these 
6  frames  is  AIRS  housekeeping  data,  the  sequencer  values  and 
status  words.  The  instrument  can  be  operated  in  any  of  three 
modes:  imaging,  spectrometer  and  photometer.  A  high  degree 
of  flexibility  exists  for  commanding  the  wavelengths  to  be 
selected  for  detectors  1  and  2  in  imaging  and  photometer 
mode.  Detectors  3  and  4  can  be  commanded  to  either  of  two 
fixed  wavelengths. 

In  imaging  mode,  a  mirror  line  scan  consisting  of  326 
pixels  takes  place  for  each  of  the  4  detectors.  Each  step  in 
the  line  scan  corresponds  to  approximately  .4  degrees. 
Although  the  instrument  actually  moves  through  336  steps,  no 
data  is  taken  in  the  first  and  last  5  steps.  Only  the  data 
from  the  326  good  pixels  is  readout  to  telemetry.  The  full 
mirror  scan  corresponds  to  plus  and  minus  67.2  degrees  from 
nadir.  At  the  end  of  each  forward  scan,  the  mirror  direction 
is  reversed  by  stepping  the  motor  at  4  times  the  forward  scan 
rate . 
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In  photometer  mode,  the  mirror  is  fixed  at  nadir.  As  in 
imaging  mode  336  pixels  are  readout  over  6  SDF  telemetry 
frames  but  valid  photometer  data  is  obtained  for  all  pixels 
except  the  first  and  last  5.  In  this  mode,  detectors  1  and  2 
are  used  for  VUV  measurements  while  detectors  3  and  4  are  for 
UV  measurements . 

Spectrometer  mode  has  the  mirror  fixed  at  nadir  while 
the  wavelength  range  is  swept.  Detector  1  covers  the 
wavelength  range  from  1189  to  2164  Angstroms  while  detector  2 
scans  the  range  from  924  to  1924  Angstroms.  Of  the  336  pixel 
values  readout  over  the  6  SDF  telemetry  frames,  valid 
spectrometer  data  is  obtained  from  pixels  6  through  331.  In 
this  mode,  detectors  3  and  4  operate  as  photometers. 

The  pixel  data  is  in  8  bit  compressed  count  form.  The 
exponent  (3  most  significant  bits)  and  mantissa  (5  least 
significant  bits)  must  be  extracted  and  used  in  the 
decompression  algorithm  in  order  to  obtain  counts  per 
accumulation  period. 

The  AIRS  data  base  was  designed  with  sufficient 
flexibility  to  be  used  in  follow  on  analyses  and  display 
routines.  The  data  base  contains  structured  records 
corresponding  to  a  full  set  of  pixel  readouts  starting  at 
pixel  1  and  ending  at  the  last  pixel. 

Decompressed  counts  are  stored  in  packed  form  to 
minimize  storage  requirements.  The  frames  are  time  tagged  to 
the  beginning  of  the  scan.  Selected  ephemeris  and  magnetic 
parameters  as  well  as  vehicle  pitch,  yaw  and  roll  are 
extracted  from  the  Time/Geometry/Magnetic  Field  records  on 
the  summary  tape  and  stored  in  each  data  base  record.  In 
addition,  the  magnetometer  data  is  stored  in  compressed  and 
packed  form  to  be  readily  available  in  the  event  it  is  needed 
for  correlation  with  the  AIRS  data.  By  storing  the 
magnetometer  data  in  compressed  counts  data  base  compaction 
is  further  enhanced. 


All  instrument  housekeeping  data  is  included  in  the  data 

base . 

Each  data  base  file  represents  data  from  one  pass .  A 
history  file  of  orbit  start  and  stop  times  along  with 
instrument  mode  is  created  using  automated  procedures.  This 
file  is  useful  in  providing  a  history  of  instrument 
operations  and  also  for  selecting  particular  passes  for 
studies . 

Since  there  is  a  great  similarity  between  the  AIRS  and 
the  AIM  sensor  flown  on  HILAT,  techniques  developed  to  handle 
the  AIM  data  were  directly  applied  to  the  AIRS.  Moreover, 
many  of  the  computer  modules  in  the  HILAT  library  of 
subroutines  have  direct  application  while  others  required 
only  minor  modification  to  have  applicability. 

The  AIRS  telemetry  consists  of  221  8--bit  words  per  SDF 
frame.  A  full  set  of  AIRS  readouts  is  acquired  over  6 
telemetry  frames.  The  first  of  the  221  words  in  each  frame 
contains  the  AIRS  frame  counter  which  cycles  from  1  to  6 
along  with  a  sequencer  A  value.  This  word  is  refered  to  as 
the  Identifier  word.  There  are  five  A  sequencers  and  one  B 
sequencer  which  are  used  to  control  the  instrument. 

In  lines  1-5,  the  telemetry  format  is  the  same:  the 
identifier  word  followed  by  220  pixel  readouts.  The  order  of 
the  pixel  data  is  as  follows  for  each  detector: 

1, 2, 3, 4, 1, 2,  .  .  . 4  .  Thus  there  are  55  readouts  for  each  of  the 
four  detectors  on  each  of  the  first  5  lines  of  telemetry. 

Line  six  begins  the  same  way  as  the  other  5  lines  but  ends  in 
a  different  way.  The  first  word  is  the  identifier  word;  the 
next  204  words  contain  pixel  data  (structured  as  in  the 
previous  5  frames);  words  206-221  contain  AIRS  housekeeping 
and  status  information.  There  are  16  AIRS  housekeeping  words 
which  are  readout  in  4  groups  of  4  words . 

The  prime  functions  of  the  A  sequencers  are: 
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LV 


Li 


AO 


Selects  one  of  16  wavelength  positions  for 
imaging  or  photometer  mode. 


A1  Selects  imaging,  photometer  or  spectrometer 

mode . 

A2  Controls  backup  functions  such  as 

increasing  power  for  spectrometer  grating 
motor . 

A3  Selects  test  modes  such  as  dark  shutter  and 

optical 

A4  Switches  power  off  to  any  of  the 

detectors 

A5  Controls  spectrometer  grating  mode 

selection  and  other  housekeeping  functions. 


The  B  sequencer  is  used  to  select  the  A  sequencer  which 
is  to  be  commanded. 

The  correlation  between  the  AIRS  frame  counter  and  the  A 
sequencer  number  is: 

Frame  Counter  A  Sequencer 

1  AO 

2  A1 

3  A2 

4  A3 

5  A4 

6  A5  . 


For  sequencer  AO,  the  sequencer  value  reflects  the 
wavelength  selected  for  detectors  1  and  2. 


0 

1216/  976 

1 

1304/1064 

2 

1356/1116 

3 

1410/1170 

4 

1456/1216 

5 

1493/1253 

6 

1544/1304 

7 

1596/1356 

8 

1625/1385 

9 

1654/1414 

10 

1670/1430 

11 

1733/1493 

12 

1750/1510 

13 

1833/1593 

14 

1910/1670 

15 

1990/1750 

Sequencer  A1  is  used  to  determine  instrument  scan  mode 


Al.  value  Mode (I=Imaging. S=Spectrometer. P=Photometer) 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


I 

S 

I 

P 

I 

s 

I 

p 

I 

s 

I 

p 

I 

s 

I 

p 


Sequencer  A3  is  used  to  determine  the  instrument  test 


0 

Normal 

1 

Dark  Shutter 

2 

Normal 

3 

Optical  Test 

4 

Normal 

5 

Extended  Dark 

6 

Normal 

7 

Optical  Test 

8 

Normal 

9 

Dark  Shutter 

10 

Normal 

11 

Optical  Test 

12 

Normal 

13 

Extended  Dark 

14 

Normal 

15 

Optical  Test 

In  addition  to  the  wavelengths  for  detectors  1  and  2 
which  are  obtained  from  sequencer  AO,  other  wavelengths  may 
be  selected.  The  selected  wavelength  can  be  determined  by  use 
of  the  A1  and  A3  sequencers. 

The  AIRS  data  base  routine  performs  the  following  prime 
funct ions : 

1.  Find  the  first  occurance  of  a  line  counter  value  of  1. 
(This  will  initiate  processing  for  the  pass.) 

2.  Test  for  telemetry  dropout  and  l's  fill  within  scans 
as  necessary.  Missing  frames  can  be  determined  by  the 
use  of  the  AIRS  line  counter  and  the  SDF  frame 
counter . 

3.  Interpolate  ephemeris/magnet ic/attitude  data  to  the 
time  at  the  beginning  of  the  line  scan. 

4.  Extract  and  decompress  pixel  counts  for  each 
detector . 

5.  Structure  data  by  detector  in  packed  word  form. 

6.  Extract  and  retain  AIRS  status  and  housekeeping  data . 
These  words  will  be  kept  in  telemetry  counts. 

7.  Extract  and  retain  data  from  the  science  magnetometer 


8.  Create  2  data  base  header  records.  The  first  contains 
a  summary  of  instrument  status  at  the  beginning  of  the 
pass.  The  second  contains  ephemeris  parameters  which 
may  be  required  in  future  AIRS  analyses. 

9.  Create  structured  AIRS  data  base  records  designed  to 
allow  for  easy  input  to  analysis  routines  yet  provide 
for  data  compression. 

The  AIRS  data  base  format  is  contained  in  the  appendix. 


4.0  Accelerometer  Analyses 

The  problem  presented  to  the  laboratory  was  to  develop 
and  implement  analytic  techniques  necessary  to  create  and 
maintain  neutral  atmospheric  density  data  bases  for  three 
satellite  missions  refered  to  as  NP5,  NP6  and  S85-1.  Analysis 
and  software  system  development  was  carried  out  in  order  to 
I  >r  ore  data  from  t.he  triaxial  accelerometer  system.';  onboard 

these  satellites.  Atmospheric  density  was  obtained  in  the 
altitude  region  from  170  km  to  240  km.  These  data  were 
merged  with  satellite  ephemeris  parameters  and  stored  in 
geophysical  history  data  bases  for  later  analysis. 

The  triaxial  accelerometer  systems  were  designed  to 
determine  atmospheric  density  by  measuring  the  satellite 
deceleration  caused  by  aerodynamic  drag.  Each  system 
consisted  of  three  accelerometer  sensors  mounted  on 
orthogonal  axes.  Thus  highly  accurate  measurements  of 
aerodynamic  accelerations  acting  on  the  satellite's  long 
track  (z),  cross  track  (x)  and  radial  direction  (y)  were 
made  . 

The  problem  consists  of  an  analysis  of  calibration  data 
for  the  accelerometers,  calculation  of  drag  coefficients  for 
the  individual  satellites,  removal  of  bias  and  noise  effects 
from  the  accelerometer  signal,  calculation  of  ephemeris 
parameters,  computation  of  atmospheric  model  values  and  their 
comparison  to  derived  values,  and  the  calculation  of  density. 

Each  acceleration  sensor  consisted  of  an 
electrostatically  suspended  proof  mass  which  was  also 
electrostatically  rebalanced  along  a  sensitive  axis  with  a 
force  equal  to  the  applied  acceleration.  Any  acceleration  of 
the  sensor  covering  caused  the  proof  mass  to  move  in  the 
opposite  direction  relative  to  the  covering.  The  rate  at 
which  pulses  of  a  given  voltage  were  required  to  restore  the 
proof  mass  to  the  center  position  was  directly  proportional 
to  the  applied  acceleration  along  the  sensitive  axis. 


For  each  axis  of  the  accelerometer,  the  output  signal 
was  generated  as  pulse  counts,  representing  the  integrated 
acceleration  over  the  sampling  period  of  about  2  seconds.  The 
calibration  data  consisted  of  tables  of  pulse  counts  per 
second  versus  acceleration  for  known  gravitational 
accelerations  applied  to  each  axis  of  the  accelerometer . 

Data  was  taken  at  different  ambient  temperatures  for  the 
accelerometer  to  determine  the  possible  temperature 
dependence  of  the  counts-to-acceleration  conversion. 


The  scale  factor  is  defined  as  the  ratio  of  counts  to 
acceleration.  Polynomial  best  fits  of  various  orders  were 
performed  for  variation  of  scale  factor  with  temperature  for 
each  of  the  various  operating  ranges. 

In  order  to  obtain  the  density  data  from  the  instrument 
output  pulse  rate,  the  following  conditions  were  required: 
the  data  was  merged  with  ephemeris  and  satellite  aspect  data; 
a  correction  was  applied  to  the  measured  pulse  rate  to 
account  for  temperature  sensitivity  and  scale  factors;  the 
pulses  due  to  air  drag  were  separated  from  those  due  to  all 
other  causes;  the  drag  acceleration  was  determined;  and  the 
density  calculated. 


For  each  axis  of  the  accelerometer,  the  measured 
acceleration  is  related  to  the  atmospheric  density  (rho)  and 
apparent  atmospheric  velocity  (V)  according  to: 


A 


rho 


2 

V  C 


4 

(V) 


2M 


(4-1) 
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where:  C  (V)  is  the  drag  coefficient  for  axis  i  (i=x,y,z), 
and  M  is  the  satellite  mass. 


Because  the  drag  coefficient  cannot  be  determined  by 
direct  measurement  for  such  low  densities  and  high 
velocities,  it  has  been  calculated  from  a  theoretical 
representation  of  the  satellite  shape.  The  computation  of 
drag  coefficients  was  dependent  on  the  ambient  temperature  of 


the  gas,  it's  mean  molecular  weight  and  the  velocity 
components  relative  to  the  accelerometer  coordinate  system. 
These  results  were  incorporated  into  subroutines  for  use  in 
the  bias  determination  and  density  calculation  stages  of  the 
processing . 

In  addition  to  the  accelerations  produced  by  atmospheric 
forces,  the  accelerometers  also  respond  to  satellite 
vibrations,  attitude  and  orbital  dynamics,  and  the 
accelerometer  output  signal  contains  noise  intrinsic  to  the 
electronics.  Because  the  output  signal  is  sampled  at  a  rate 
of  about  0.5  Hz,  the  extraneous  signals  can  make  a 
contribution  to  the  slowly  varying  atmospheric  signal  through 
aliasing,  although  some  attenuation  of  aliased  signal  does 
occur  because  of  the  2-second  integration  time  per  sample. 
High-frequency  noise  is  further  reduced  by  an  electronic 
filter  associated  with  the  accelerometer  electronics.  This 
is  done  prior  to  sampling  so  that  contributions  at  the 
aliased  frequency  are  also  eliminated. 

The  attitude  adjustments  and  vibrations  were  treated  by 
low  pass  digital  filtering  of  the  various  accelerometer 
signals.  The  parameters  chosen  for  the  low-pass  filters  for 
the  individual  axes  were: 


x-axis 

y-axis 

z-axis 


Pass  Bandwidth  Transition  Bandwidth 


0.05  Hz 
0.05  Hz 
0.05  Hz 


0.0167  Hz 
0.0167  Hz 
0.0167  Hz 


In  the  case  of  the  NP5  vehicle  an  additional  data  base 
was  created  using  a  low-pass  filter  with  the  following 
characteristics : 

Pass  Bandwidth  Transition  Bandwidth 

x-axis  0.005  Hz 

y-axis  0.005  Hz 

z-axis  0.010  Hz 


0.005  Hz 
0.005  Hz 
0.010  Hz 


Extraordinary  noise  contributions  could  not  be  treated 
by  digital  filtering,  and  were  of  such  a  magnitude  as  to 
overwhelm  the  aerodynamic  signal.  Data  segments  affected  by 
these  were  simply  skipped  when  density  values  were 
calculated . 

One  characteristic  of  the  accelerometers  was  a  non-zero 
response  to  zero  aerodynamic  acceleration  while  in  orbit. 

The  bias  values  are  different  for  each  of  the  three 
accelerometer  axes  and  are  not  constant  as  a  function  of 
time.  Thus,  a  method  of  removing  these  bias  values  from 
acceleration  signal  was  needed. 

For  each  pass,  atmospheric  drag  at  satellite  apogee 
contributed  to  the  accelerometer  signal.  This  contribution 
was  estimated  by  using  a  model  atmospheric  density  (Jacchia 
1971)  and  equation  (4-1)  to  compute  expected  drag 
accelerations  for  each  accelerometer  axis  in  the  vicinity  of: 
apogee.  For  each  axis,  the  bias  was  calculated  as  the 
difference  between  the  measured  signal  from  the  accelerometer 
and  the  corresponding  expected  drag  acceleration.  Thusly, 
averaged  biases  were  taken  over  various  time  periods  centered 
on  apogee . 

The  resulting  bias  values  showed  some  scatter  about  a 
much  slower  trend  over  each  vehicle's  lifetime.  This  scatter 
was  attributed  to  the  estimations  of  the  atmosphere  model  and 
noise  in  the  accelerometer  signal,  after  filtering.  Thus,  a 
smooth  fit  to  the  data  was  used  to  represent  the  bias  value. 
For  each  vehicle  and  accelerometer  axis,  polynomial  fits  wen- 
performed  over  vehicle  lifetime. 

Accelerometer  data  was  merged  with  orbital  ephemeris 
data  based  on  corresponding  time  values,  although  orbital 
data  was  interpolated  from  values  spaced  at  60  second 
intervals  to  the  roughly  2  second  interval  of  the 
accelerometer  data.  A  quadratic  interpolation  process  was 
used  for  this  evaluation  for  the  fundamental  ephemeris 
parameters  (e.g.,  latitude,  longitude,  altitude)  and  ah;-  f  t 


some  derived  parameters  (e.g.,  heading  direction, 
velocities).  Other  derived  parameters  were  computed  on  a 
point-  by-point  basis  from  more  fundamental  parameters  (e.g., 
magnetic  local  time,  local  time) . 

The  following  parameters  were  added  to  the  accelerometer 
data  base  using  the  supplied  ephemeris  values: 

a)  altitude 

b)  latitude 

c)  longitude  (positive  east) 

d)  geocentric  velocity  (km/sec) 

e)  local  time 

The  remaining  parameters  were  derived  from  the  ephemeris 
values  : 

a)  velocity  with  respect  to  atmosphere  (km/sec) 

b)  angle  between  geocentric  and  atmospheric  velocities 

c)  solar  declination 

d)  x-direction  relative  to  North 

e)  magnetic  local  time  . 

Density  values  were  derived  from  the  measured 
accelerations  using  the  following  equation, 

2  M  Az 

rho  - -  ,  where 

2  + 

V  Cz (V) 

Az  is  the  aerodynamic  acceleration  for  axis  z  (long 
track) , 

M  is  the  satellite  mass, 

V  is  the  total  atmospheric  velocity  relative  to  the 
satellite, 

V  is  the  satellite  velocity, 

Cz  (V)  is  the  drag  coefficient  for  axis  z,  and 
rho  is  the  atmospheric  density. 

The  total  atmospheric  velocity  relative  to  the  satellite 
is  given  by 

■*  ■*  +  -» 

V  =  -Vg+Va+Vw,  where 

Vg  =  geocentric  velocity  of  the  satellite, 


Va  =  rotational  contribution  of  atmospheric  velocity,  and 

Vw  =  wind  contribution  of  atmospheric  velocity. 

For  this  study,  the  wind  contribution  being  small  in 
relation  to  the  known  contribution  -Vg+Va,  was  ignored. 

A  set  of  existing  computer  routines  were  modified  and 
adapted  to  process  each  system,  uniquely.  In  the  case  of  NFS 
and  NP 6,  digital  tapes,  refered  to  as  DMA  tapes,  containing 
raw  and  filtered  accelerometer  data  were  received  at  AFGL . 

In  the  case  of  S85-1,  LPARL  tapes  containing  raw 
accelerometer  data  were  provided.  The  data  was  processed  on  a 
daily  basis  to  the  final  density  stage.  Associated  displays, 
listings  and  data  bases  were  created  as  noted  below. 

For  each  day  in  question,  the  filtered  and  raw  data  was 
reformatted  onto  a  magnetic  tape.  In  the  case  of  S85-1,  the 
raw  data  was  filtered  initially. 

Time  history  displays  were  produced  which  exhibit  the 
data  for  each  of  the  3  accelerometer  axes,  along  with 
temperature  on  a  daily  basis.  The  data  was  then  merged  with 
appropriate  ephemeris  parameters  and  stored  on  a  magnetic 
tape.  A  decimated  merged  data  base  was  stored  on  permanent 
file  for  further  plotting.  The  merged  data  along  with 
appropriate  Kp,  Ap  and  flux  values  was  accessed  in  the 
generation  of  bias  values  for  the  3  axes  on  an  orbit-by-orbit 
basis  and  a  daily  average  calculation.  The  bias  calculations 
were  based  on  fitting  the  filtered  data  to  model  drag 
components  over  500  seconds  of  apogee  data  for  each  orbit . 

In  the  case  of  S85-1,  180  seconds  of  apogee  data  were  used  in 
the  averaging  process.  Paper  plots  and  listings  were 
generated  for  each  day  in  question.  Also,  the  decimated 
merged  data  file  was  accessed  in  the  generation  of  paper 
plots  exhibiting  filtered  outputs  for  each  of  the  3  axes  te: 
an  orbit-by-orbit  basis.  At  this  point,  the  bias  filer, 
(orbit-by-orbit  and  daily  average  files)  augmented  by  the  day 
being  processed  were  updated.  The  files  were  sorted  day-by- 


day,  EOR ' s  removed  and  duplicate  days  eliminated  prior  to 
usage  by  the  density  program. 


m 


On  the  basis  of  'quiet'  vs.  all  days,  plots  exhibiting 
bias  values  as  a  function  of  time  were  generated  resulting  in 
a  bias  fit  technique  being  used  for  density  calculations. 
Based  on  a  study  of  the  raw  and  filtered  data  plots,  skip 
times,  necessary  to  the  density  program  run,  were  determined 
for  each  day.  A  density  data  base  was  created  on  tape  and  on 
permanent  file  with  the  JACCHIA  program.  Paper  plots  and 
listings  were  generated  from  the  density  file  exhibiting 
density  and  normalized  density  and  selected  model  and 
ephemeris  parameters  as  functions  of  time  on  an  orbit-by- 
orbit  basis.  Microfiche  displays  of  density  as  a  function  of 
altitude  were  generated.  In  the  case  of  S85-1,  microfiche 
displays  were  not  required.  Finally,  follow  on  density  plots 
(normalized  with  respect  to  the  MSIS  model)  were  generated 
for  a  wave  structure  study. 

Raw  and  filtered  digital  data  were  provided  for  156  days 
of  NP5  and  160  days  of  NP 6 .  A  total  of  150  days  of  NP 5  and 
150  days  of  NP 6  were  processed  to  the  final  density  stage. 

The  remaining  6  days  of  NP 5  and  10  days  of  NP 6  were 
unprocessable . 

For  S85-1,  78  days  of  raw  data  were  provided  for 
processing.  Each  day  contained  approximately  ten  28  minute 
passes  centered  on  maximum  north  latitude,  six  20  minute 
passes  centered  on  perigee  and  eight  6  minute  passes  centers  : 
on  apogee  for  bias  determination.  For  each  data  base  tape, 
accelerometer  words,  clocks  and  marker  bits  were  stripped 
off,  and  the  data  was  time  correlated  and  reconstructed  in 
proper  time  sequence.  The  accelerometer  words  were  then 
calibrated  to  mg's  and  temperature  and  data  base  files  were 
created.  For  this  purpose  calibrations  and  clock  correlation 
factors  were  utilized.  Analyses  and  associated  computer 
routines  were  developed  and  successfully  utilized  to  convert 
^  he  time  correlated  raw  data  to  filtered  data.  Raw  and 


filtered  data  tapes  were  thereby  generated  in  the  same  format 
as  the  reformatted  DMA  tapes  for  the  78  days  of  interest .  A 
total  of  74  of  the  78  days  were  subsequently  processed  to  the 
final  density  stage.  The  remaining  4  days  were  unprocessable. 

The  density  data  bases,  noted  above,  were  generated  at 
high  data  rate  outputs  (a  point  every  4  seconds) .  These  data 
bases  provided  the  opportunity  to  store  data  from  many  orbits 
into  common  data  banks  which  were  stuctured  to  allow  global 
studies  to  be  performed. 

Tape  formats  associated  with  raw/filtered  accelerometer 
data,  merged  ephemeris  data  and  model/density  data  are 
contained  in  the  appendix. 
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5.0  Shuttle  Analyses 


Analysis  systems  have  been  developed  and  implemented  for 
space  shuttle  missions.  Data  bases  consisting  of  experiment 
data  along  with  shuttle  related  parameters  have  been 
successfully  created.  Systematic  techniques  analagous  to 
those  described  in  preceding  sections  were  employed  for  the 
experiment  data.  New  techniques  were  necessary  for  the  data 
products  received  at  AFGL  from  NASA.  This  chapter  will 
discuss  the  NASA  products. 

For  early  missions  two  sets  of  tapes  were  typically 
received.  The  preliminary  data  (on  tapes  refered  to  as 
CP011F)  was  received  shortly  after  landing.  The  CP011F  tapes 
contained  unprocessed  orbiter  data  as  it  was  downlinked  over 
several  asynchronous  PCM  systems.  Time  tagging  was 
accomplished  by  using  appropriate  delta  time  values  contained 
on  the  records.  The  final  tapes  (CBET04)  contained  smoothed 
and  processed  data  but  these  tapes  begin  to  arrive 
approximately  12  weeks  after  the  mission.  On  these  two  tape 
types,  the  orbit,  attitude  and  ancillary  data  are  blocked 
together  in  large  records. 

For  recent  missions,  the  NASA  data  has  been  received  on 
Shuttle  POCC  Interface  Facility  (SPIF)  tapes.  These  tapes  are 
similar  to  a  combination  of  CP011F  and  CBET04  tapes  but  the 
orbit,  attitude  and  ancillary  data  is  interlaced  on  each 
file.  Moreover,  the  tapes  were  generated  on  a  VAX  computer 
and  the  standard  byte  reversal  techniques  had  to  be  applied 
before  the  informational  parameters  could  be  extracted. 
Analytic  techniques  developed  for  the  CP011F  and  CBET04  tapes 
had  direct  application  for  the  SPIF  tapes. 

Routines  and  techniques  were  developed  for  the 
extraction  of  parameters  from  these  NASA  products.  This, 
however,  is  frequently  only  the  first  phase  of  the  effort 
since  additional  analyses  must  be  performed  in  order  to 
provide  the  experimenter  with  the  parameters  which  are 
requ i red . 
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With  respect  to  the  ephemeris  data,  the  preliminary 
Johnson  Space  Center  (JSC)  CP011F  tape  contains  two  sets  of 
vehicle  position  and  velocity  vectors  (referred  to  as  state 
vectors)  in  an  Earth  Centered  Inertial  (ECI)  true  of  date 
coordinate  system.  This  coordinate  system  (XYZ)  is  right 
handed,  where  X  is  the  component  in  the  direction  of  the 
Vernal  Equinox;  Z  is  in  the  direction  of  the  North  pole;  and 
Y  =  ZxX.  One  set  of  provided  vectors  consists  of  the 
filtered  components  outputted  by  the  telemetry,  and  the 
others  are  the  current  state  vectors  of  the  spacecraft.  The 
CBET04  JSC  tape  contains  processed  information.  The 
spacecraft's  state  vector  is  provided  in  many  coordinate 
systems  such  as  mean  of  50  ECI,  true  of  date  ECI,  and  Aries 
true  of  date  polar. 

Existing  ephemeris  computation  routines  were  then  used 
to  produce  a  standardized  ephemeris  file.  This  ephemeris  file 
is  created  at  a  sufficiently  high  sample  rate  so  that  it  can 
be  interpolated  upon  whenever  necessary  without  loss  of  data 
structure . 

The  direct  correlation  between  ephemeris  parameters  and 
generated  geomagnetic  and  solar  models  is  well  known.  Since 
experimenters  often  require  magnetic  pitch  angles  to  properly 
interpret  their  sensor  data,  various  geomagnetic  parameters 
have  to  be  determined.  Thus,  using  a  standardly  formatted 
ephemeris  file  as  input,  the  model  geomagnetic  field 
parameters  can  be  determined  using  the  appropriate  set  of 
IGRF  eighth  degree  Gaussian  coefficients.  The  proper  set 
depends  upon  the  year  of  the  flight.  For  each  time  related 
value,  the  cartesian  components  of  the  Earth's  magnetic  field 
(North,  East,  and  vertical) ,  the  magnetic  declination  and 
inclination  (dip)  angles,  can  be  determined.  The  spatial 
parameters  commonly  required  include  the  total  magnetic  field 
intensity  and  the  L  shell  (in  Earth  radii) .  Other 
geomagnetic  parameters  which  can  be  computed  are  geomagnetic 


local  time,  geomagnetic  latitude  and  longitude,  and  corrected 
geomagnetic  latitude,  longitude,  and  local  time. 

For  certain  sensors  meaningful  interpretation  of  the 
instrument's  output  can  be  accomplished  only  when  the 
orientation  of  the  sensor  is  known.  The  orientation  of  the 
sensing  axis  of  an  experiment  with  respect  to  a  fixed 
coordinate  system  in  space  is  referred  to  as  its  attitude. 

The  NASA  CP011F  tape  provides  vehicle  attitude 
information  in  the  form  of  quaternions  which  relate  the 
vehicle's  axes  in  the  M50  coordinate  system  to  the  local 
vertical,  local  horizontal  (LVLH)  coordinate  system.  Vehicle 
pitch,  yaw,  and  roll  in  the  M50  coordinate  system  are  also 
provided . 

The  CBET04  tape  contains  processed  and  filtered  results 
at  a  fixed  data  rate.  This  tape  contains  attitude  of  the 
vehicle  axes  in  quaternion  or  matrix  form.  This  data  is 
given  in  the  M50  coordinate  system.  The  attitude  of  the 
shuttle  axes  is  given  in  the  form  of  3x3  matrices  or 
quaternions.  Euler's  theorem  states  that  any  finite  rotation 
of  a  rigid  body  can  be  expressed  as  a  rotation  through  some 
angle  about  a  fixed  axis,  namely  a  quaternion.  Thus,  any 
three  axis  Euler  rotation  3x3  matrix  can  be  described  as  a 
rotation  by  some  angle  about  some  fixed  axis. 


A  three  angle  Euler  rotation,  which  relates  system  A  to 
system  B,  can  be  represented  by  a  3x3  rotation  matrix  R 
comprised  of  elements  R^j  such  that 
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The  quaternion  is  a  vector  with  four  components  which 
relates  system  A  to  system  B.  It  is  denoted 

QBA  =  (qi,q2,q3,q4) 
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The  rotation  angle  and  the  single  Euler  axis  that 
defines  the  quaternion  relationship  between  the  two  systems 
can  be  determined  from  the  3x3  matrix  components.  The 
rotation  angle  is  given  by 

cosa  ~  1/2  +R22  +*33  —  1.)  , 

and  the  single  Euler  axis 

E  =  (EX,  EY,  EZ ) 

where 

EX  —  (R23-R32 )  / 2  sina 

EY  =  (R31-R13)  / 2  sina 

EZ  =  (Ri2~R2i)  /2  sina 

Thus  QBA  =  (qi, q2» q3/ q4)  is  given  by 

QBA  =  (cos (a/2) , EX  sin(a/2),  EY  sin(a/2),  EZ  sin (a/2))  . 

Similarly,  the  3x3  rotation  matrix  can  be  determined 
from  the  quaternion.  In  this  case,  the  transpose  must  be 
taken  in  order  to  obtain  the  proper  matrix.  The  3x3  matrix 
determined  from  QBA  =  (qif  q2r  ^3^  ^4)  is 


*11 

<3131  + 

q2*32  -  q3q3  ~  q4q4 

*12 

= 

2  (q2c?3 

-qiq4> 

*13 

= 

2  (q2c?4 

+  qiq3> 

*21 

= 

2  (q2q3 

+  qiq4> 

R22 
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r33 
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qiqi  - 

q2q2  ~  q3q3  +  q4q4 
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The  most  basic  coordinate  system  used  is  the  body  axis 
coordinate  system.  In  this  system  X  is  parallel  to  the 
orbiter  structural  axis  positive  toward  the  nose,  Y  is 
parallel  to  the  right  wing,  and  Z  is  positive  out  of  the 
bottom  of  the  orbiter.  The  Euler  rotation  associated  with 
this  system  is  a  pitch,  yaw,  roll  sequence. 

The  vehicle  body  system  is  frequently  required  to  be 
expressed  in  the  Local  Vertical  Local  Horizontal  (LVLH) 
coordinate  system  (or  landing  reference  frame) .  This  is  a 
local  orbital  coordinate  system  where  the  Z  axis  goes  from 
the  vehicle  to  center  of  the  Earth,  Y  axis  is  normal  to  the 
orbital  plane  and  opposite  sense  of  the  vehicle's  momentum 
vector,  and  X  =  Y  x  Z. 

For  attitude  given  in  the  Mean  1950  (M50)  coordinate 
system,  X  is  directed  towards  the  mean  vernal  equinox  of 
1950,  the  Z  is  in  the  direction  of  the  mean  1950  North 
celestial  pole,  and  Y  completes  the  right  hand  coordinate 
system.  Another  frequently  required  attitude  reference  frame 
is  the  true  of  date  ECI  coordinate  system.  In  this 
coordinate  system  X  is  the  direction  of  the  actual  vernal 
equinox,  Z  is  in  the  direction  of  the  true  North  celestial 
pole,  and  Y  completes  the  right  hand  system. 

Since  different  parameters  can  be  expressed  in  various 
coordinate  systems,  a  series  of  matrix  transformations  has 
been  developed  to  allow  for  mappings  between  systems.  These 
routines  have  been  developed  in  modular  form  to  keep  within 
the  overall  SSDAS  philosophy. 

Transformation  matrices  relating  the  M50  reference  frame 
to  the  true  of  date,  the  true  of  date  to  the  LVLH,  and  the 
true  of  date  to  the  local  cartesian  system  exist. 

The  (XL,YL,ZL)  axes  of  the  LVLH  are  determined  from  the 

true  of  date  state  vectors.  R  is  the  position  of  the  shuttle 
from  the  center  of  the  Earth  in  ECI  system  and  it  is  given  by 

R  =  iRl  +  jR2  +  kR3 
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Let  the  absolute  value  of  a  vector,  R,  be  denoted  by  [R] , 
then 


[R]  =  VRl2  +  R22  +  R32 

Then  expressing  R  as  a  unit  vector, 

R  =R/[R]  =  iRl/tR]  +  jR2/[R]  +  kR3/ [R]  . 

The  vehicle's  velocity  vector  is  given  by 

V  =  iVl  +  jV2  +  kV3 
where  [V]  is  defined  similarly. 

In  the  LVLH  coordinate  system 

ZL  =  -R/[R] 

Yl  =  -  (RxV)  /  [  (RxV)  ] 

X  ZL 

Then  (XL,YL,ZL)  in  matrix  form  is  given  by 


Taking  the  transpose 


Thus,  by  substituting,  the  vehicle  axes  are  now 
represented  in  the  LVLH  coordinate  system  by 


Similar  transformation  matrices  have  been  developed  for 
the  determination  of  shuttle  pitch,  yaw  and  roll. 

The  attitude  of  the  LOS  of  the  vehicles'  axes  is 
normally  determined  in  all  the  above  coordinate  systems  for 
the  CP011F  and  CBET04  data.  Additionally,  the  LOS  of  probes 
aboard  the  shuttle  are  normally  computed. 

Attack  angles  and  tangent  height  calculations  are 
frequently  required  by  experimenters  for  proper  data 
interpretation . 


Attack  angles  are  defined  as  the  angles  between  the  LOS 
of  a  probe  or  axis  and  a  given  vector.  The  more  frequently 
required  attack  angles  involve  the  angle  between  a  sensor  LOS 
and  the  magnetic  field,  velocity  vector  or  the  sun  vector. 

Tangent  height  is  the  minimum  altitude  above  the  surface 

of  the  Earth  which  an  unrefracted  light  ray  coming  from 
behind  the  Earth  through  the  CO2  layer  to  the  spacecraft . 

Thus,  tangent  height  is  apparent  height  (at  the  horizon)  from 
which  radiation  is  coming.  Tangent  height  calculations  assume 
an  oblate  spheroid  earth  model  and  is  a  geometric  solution 
rather  than  the  method  of  closest  approach.  These 
calculations  also  require  attitude  of  the  LOS  in  the  local 
cartesian  system,  (elevation  and  azimuth  angles),  the 
vehicle's  positional  data  of  latitude,  altitude  and  distance 
from  the  center  of  Earth  to  vehicle,  and  an  equation  to 
calculate  the  radius  of  the  Earth  as  a  function  of  the 
geodetic  latitude. 

A  final  irreducible  data  base  consisting  of  attitudinal, 
positional,  and  geomagnetic  parameters  was  created  for  each 
mission.  This  data  base  was  designed  to  include  all  the 
parameters  necessary  for  the  satisfaction  of  experimenter 
requirements.  Moreover,  the  data  base  is  constructed  with  an 
even  time  spacing  to  facilitate  interpolation . 
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With  respect  to  the  ancillary  data,  techniques  and 
computer  routines  have  been  developed  for  the  determination 
of  thruster  firings,  occurances  of  water  dumps,  OMS  burns. 
These  parameters  as  well  as  others  such  as  fuel  cell  purges 
and  relief  of  the  cryogenic  supplies  are  considered  to  be 
contaminants  since  they  can  drastically  effect  sensor 
outputs.  Typically,  a  full  mission  data  base  of  thruster 
firings  and  water  dumps  was  created. 

Onboard  the  orbiter  there  are  38  main  and  6  Vernier 
thrusters  which  are  included  in  the  Reaction  Control  System 
(RCS) .  Twelve  main  jets  and  2  Vernier  thrusters  are  located 
forward,  and  24  main  jets  (2  systems  of  12  pairs)  and  4 
Vernier  thrusters  are  located  in  the  aft  section.  Knowledge 
of  precisely  which  thruster  (or  thrusters)  were  fired  at  any 
instant  of  time  can  be  essential  to  data  interpretation. 

There  are  two  water  storage  systems  onboard  the  vehicle. 
One  is  a  supply  water  system  (which  consists  of  six  tanks) 
and  the  other  is  a  waste  water  system  (one  tank) .  Only  one 
tank  can  be  dumped  at  a  time.  The  waste  water  tank  is  full 
at  launch  and  is  dumped  in  orbit,  and  then  the  tank's  water 
level  is  restored  to  its  condition  at  launch.  Bi-levels 
within  the  spacecraft's  telemetry  indicate  the  occurrance  and 
type  of  water  dump  in  progress.  Analog  monitors  onboard  the 
spacecraft  indicate  the  percentage  of  water  in  each  tank . 
Knowing  the  capacity  of  each  tank  allows  the  calculation  of 
the  amount  of  water  being  dumped.  Thus,  for  any  water  dump 
the  following  information  can  be  determined;  the  type  of 
water  dump  (waste  or  supply) ;  the  start  and  end  time  of  each 
dump;  the  duration  and  quantity  of  water  discharged  at  any 
instant  while  the  dump  is  in  progress;  and  the  total  amount 
of  water  expelled  during  the  dump. 

Data  bases  structured  for  use  on  the  AFGL  Control  Data 
Cyber  computer  have  been  archived. 
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HILAT  TAPE  LOG  SAMPLE 
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HILAT 


SUMMARY  TAPE  FORMAT 

TIME/GEOMETRY/MAGNETIC  FIELD  SUMMARY  DATA 
330  16-bit  word  record  (note  1) 

(one  per  15-second  summary  block) 

Number  of 


Parameter  16-bit  Words 

-  Universal  Year  1 

time  Day  1 

Hour  1 

Minute  1 

Second  1 

Millisecond  1 

-  Station  Geodetic  latitude — degree  +N  2 

Longitude--degree  +E  2 

Geodetic  altitude — km  2 

NORAD  ephemeris  elements,  line  1  35 

NORAD  ephemeris  elements,  line  2  35 

Azimuth  angle  to  satellite — degree  2 

Elevation  angle  to  satellite — degree  2 

Range  to  satellite — km  2 

Local  earth  radius--km  2 

(spares)  48 

-  Satellite  Geodetic  latitude — degree  +n  2 

Longitude--degree  +E  2 

Geodetic  altitude — km  2 

Geocentric  position  X--km  (note  2)  2 

Geocentric  position  Y--km  (note  2)  2 

Geocentric  position  Z--km  (note  2)  2 

Geocentric  velocity  X--km/s  (note  2)  2 

Geocentric  velocity  Y--km/s  (note  2)  2 

Geocentric  velocity  Z--km/s  (note  2)  2 

Orbital  attitude  pitch--degree  (note  3)  2 

Orbital  attitude  yaw — degree  (note  3)  2 

Orbital  attitude  roll-degree  (note  3)  2 

Invariant  magnetic  latitude, 

APL  eccentric  dipole--degree  2 

Magnetic  local  time, 

APL  eccentric  dipole — seconds  2 

Magnetic  dip,  IGRF80--degree  2 

Magnetic  declination,  IGRF80--degree  2 

Magnetic  strength  north,  IGRF80--mG  2 

Magnetic  strength  east,  IGRF80-mG  2 

Magnetic  strength  vertical,  IGRF80--mG  2 
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hilat 

SPMMARY  TAPE  FORMAT  (Cont.) 

Number  of 

Parameter  16-bit  Words 


Satellite  Solar  right  ascension--degree  2 

Solar  decliantion — degree  2 

Solar  zenith  angle--degree  2 

Sun-shade  angle — degree  2 

Invariant  magnetic  latitude, 

Gustafsson  model — degree  2 

Magnetic  longitude, 

Gustafsson  model — degree  2 

(spares)  4 

Magnetic  field  Geodetic  latitude  (350  km) --degree  +N  2 

line  through  Longitude  (350  km) --degree  +E  2 

satellite  Geodetic  altitude  (350  km) — km  2 

(IGRF80)  Solar  zenith  angle  (350  km) --degree  2 

Sun-shade  angle  (350  km) --degree  2 

Geodetic  latitude  (100  km) — degree  +N  2 

Longitude  (100  km) --degree  +E  2 

Geodetic  altitude  (100  km) --km  2 

Solar  zenith  angle  (100  km) --degree  2 

Sun-shade  angle  (100  km) — degree  2 

Geodetic  latitude  (surface) --degree  +N  2 
Longitude  (surface) — degree  +E  2 

Geodetic  altitude  (surface) --km  2 

Solar  zenith  angle  (surface) — degree  2 

Sun-shade  angle  (surface) -degree  2 

(spares)  4 

F-region  Geodetic  latitude — degree  +N  2 

penetration  Longitude--degree  +E  2 

(350  km)  Geodetic  altitude — km  2 

Invariant  magnetic  latitude, 

APL  eccentric  dipole — degree  2 

Magnetic  local  time, 

APL  eccentric  dipole--seconds  2 

Magnetic  dip,  IGRF80--degree  2 

Magnetic  declination,  IGRF80--degree  2 

Magnetic  strength  north,  IGRF80--mG  2 

Magnetic  strength  east,  IGRF80 — mG  2 

Magnetic  strength  vertical,  IGRF80--mG  2 

Propagation  zenith  angle — degree  2 

Propagation  magnetic  azimuth--degree  2 

Off-magnetic  field  propagation  angle, 

IGRF80— degree  2 

Off-magnetic  meridian  propagation  angle, 
IGRF80--degree  2 

Off-magnetic  L-shell  propagation  angle, 

IGRF80--degree  2 


HILAX 

SUMMARY  TAPE  FORMAT  (Cont.) 

Number  of 

Parameter  16-bit  Words 

F-region  Reduced  propagation  range--km  2 

penetration  Magnetic  penetration  velocity 

(350  km)  X — km/s  (note  4)  2 

Magnetic  penetration  velocity 

Y — km/s  (note  4)  2 

Magnetic  penetration  velocity 

Z — km/s  (note  4)  2 

(spares)  6 

F-region  Geodetic  latitude--degree  +N  2 

penetration  Longitude--degree  +E  2 

(100  km)  Geodetic  altitude — km  2 

Invariant  magnetic  latitude, 

APL  eccentric  dipole — degree  2 

Magnetic  local  time, 

APL  eccentric  dipole — seconds  2 

Magnetic  dip,  IGRF80 — degree  2 
Magnetic  declination,  IGRF80 — degree  2 

Magnetic  strength  north,  IGRF80--mG  2 

Magnetic  strength  east,  IGRF80 — mG  2 

Magnetic  strength  vertical,  IGRF80--mG  2 

Propagation  zenith  angle — degree  2 

Propagation  magnetic  azimuth — degree  2 

Off-magnetic  field  propagation 

angle,  IGRF80--degree  2 

Off-magnetic  meridian  propagation 

angle,  IGRF80 — degree  2 

Off-magnetic  L-shell  propagation 

angle,  IGRF80--degree  2 

Reduced  propagation  range--km  2 

Magnetic  penetration  velocity 

X — km/s  (note  4)  2 

Magnetic  penetration  velocity 

Y — km/s  (note  4)  2 

Magnetic  penetration  velocity 

Z — km/s  (note  4)  2 

(spares)  6 

Magnetic  field  Geodetic  latitude  (100  km) --degree  +N  2 

line  thourgh  Longitude  (100  km) --degre  +E  2 

F-region  Geodetic  altitude  (100  km) --km  2 

penetration  (spares)  2 
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hilat 


note  1  -one  word  parameters  are  16-bit  integers,  two  word 
parameters  are  HP  32-bit  real  words,  and  the  NORAD 
elements  are  ASCII  coded  FORTRAN  F10.5  format 

note  2  -  geocentric  coordinates  are  right-hand  earth-centered 

inertial  with  X  through  the  line  of  Aries  and  Z  through 
the  north  pole 

note  3  -attitude  procesing/error  flags  are  as  follows: 

-999  no  attitude  solution  attempted 

-998  telemetry  frame  error 

-997  all  sensors  dark 

-996  solar  sensor  word  unpack  error 

-995  model-data  magnetic  field  magnitude  error 

-994  model-data  magnetic  field/sunline  angle  error 

-993  orthogonalization  error 

-992  orthogonalization  error 

-991  orthogonalization  error 

-990  miscellaneous  error 

note  4  -penetration  magnetic  coordinates  are  right-hand,  X 
magnetic  north  and  Z  nadir 


hilat 


-  Unpacked  TM  (sequence  repeated  31  times)  7812 


Frame  count  (note  3)  (2) 

Overrange  bit  and 

sync  status  (note  4)  (1) 

SDF  timer  (1) 

NSD/SAD  (1) 

Beacon  housekeeping  (1) 

Magnetometer  (31) 

Particle  detector  (97) 

AIM  (62) 

Drift  meter  (56) 


-  Ancillary  data  (spares) 


8 


rote  1  -  all  parameters  are  16-bit  integer  words,  except  as  noted 

note  2  -  times  are  with  respect  to  the  time/geometry/magnetic  fiel 
point  in  Table  3 

r.  'te  3-  frame  count  is  in  32-bit  integer  format 
note  4-  overrange  bit  and  sync  status  is  as  follows: 

0  frame  sync,  no  overrange 

1  frame  sync,  overrange 

2  frame  sync  error,  no  overrange 

3  frame  sync  error,  overrange 

4  no  frame  sync,  no  overrange 

5  no  frame  sync,  overrange 

6  replacement  for  missing  frame 


HUAI 


SCINTILLATION  SUMMARY  DATA 

990  16-bit  word  record  (note  1) 
(one  per  15-second  summary  block) 


Number  of 


Parameter  16-bit  words 

-  Span  Start  time  of  data — second  (note  2)  2 

End  time  of  data--second  (note  2)  2 

(spares)  6 

-  137  MHz  Mean  signal  level--db  2 

Intensity  scintillation  index  S4  2 

Intensity  fade  period--second  2 

Intensity  decorrelation  time--second  2 
Standard  deviation  of  phase--rad  2 

Intensity  power  spectral  density 
function  decimated  in  log 
f requency--db  S4  units  sq/Hz  100 

Phase  power  spectral  density 
function  decimated  in  log 
f requency--db  rad  sq/Hz  100 

Frequencies  of  above  power 

spectral  density  samples--Hz  100 

-  413  MHz  Mean  signal  level — db  2 

antenna  1  Intensity  scintillation  index  S4  2 

Intensity  fade  per iod--second  2 

Intensity  decorrelation  time-second  2 

Standard  deviation  of  phase--rad  2 

Intensity  power  spectral  density 
function  decimated  in  log 
f requency--db  S4  units  sq/Hz  100 

Phase  power  spectral  density 
function  decimated  in  log 
f requency--db  rad  sq/Hz  100 

Frequencies  of  above  power 

spectral  density  samples — Hz  100 

-  1239  MHz  Mean  signal  level--db  2 

Intensity  scintillation  index  S4  2 

Intensity  fade  period--second  2 

Intensity  decorrelation  time--second  2 
Intensity  power  spectral  density 
function  decimated  in  log 
f requency--db  S4  unts  sq/Hz  100 

Frequencies  of  above  power 

spectral  density  samples--Hz  100 


EILAT 


SCINTILLATION  SUMMARY  DATA  (Cont.) 


30  MHz 


Parameter 

Mean  signal  level--db 
Intensity  scintillation  index  S4 
Intensity  fade  period--second 
(spare) 

Standard  deviation  of  phase--rad 
(spares) 


Number  of 
16-bit  words 

2 

2 

2 

2 

2 

10 


113  MHz 
antenna  2 


413  MHz 
antenna  3 


-  435  MH: 


:;-:c 


Irregularity 
anisotropy 
and  drift 


Mean  signal  level--db  2 

Intensity  scintillation  index  S4  2 

Intensity  fade  period--second  2 

Intensity  decorrelation  time--second  2 
Standard  deviation  of  phase--rad  2 

(spares)  10 

Mean  signal  level — db  2 

Intensity  scintillation  index  S4  2 

Intensity  fade  period--second  2 

Intensity  decorrelation  time — second  2 
Standard  deviation  of  phase--rad  2 

(spares)  10 

Mean  signal  level — db  2 

Intensity  scintillation  index  S4  2 

Intensity  fade  period--second  2 

(spare)  2 

Standard  deviation  of  phase--rad  2 

(spares)  10 

Total  electron  content  along 
propagation  raypath,  three 
equispaced  samples--el/m  sq  6 

Axial  ratio  of  receiver  plane 

correlation  surface  2 

Geographic  azimuth  of  receiver 

plane  correlation  surf ace--degree  2 

Receiver  plane  correlation  surface 

velocity,  magnetic  north--km/s  2 

Receiver  plane  correlation  surface 

velocity,  magnetic  east--km/s  2 
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HUAI 


SCINTILLATION  SUMMARY  DATA  (Cont.) 


Parameter 

Number  of 
16-bit  words 

-  Angular 

Standard  deviation  of  phase 

deviation 

difference,  antennas  l,2--rad 

2 

Standard  deviation  of  phase 

difference,  antennas  l,3--rad 

2 

Standard  deviation  of  phase 

difference,  antennas  2,3 — rad 

2 

-  Ancillary  data 

(spares) 

52 

note  1  -  all  parameters  are  HP  32-bit  real  words;  items  not 
calculated  are  dummy  filled  with  -999 


note  2-  times  are  with  respect  to  the  time/geometry/magnetic  field 
point  and  will  always  be  -15.  and  +15.  seconds 


HILAT 

SCIENCE  SUMMARY  DATA 

1875  16-bit  word  record  (note  1) 
(one  per  15-second  summary  block) 


Number  of 

Parameter  16-bit  Words 


Span  Start  time  of  data--second  (note  2)  2 

End  time  of  data — second  (note  2)  2 

(spares)  6 

Vector  X-axis  data — counts  30 

Magnetometer  Y-axis  data — counts  30 

Z-axis  data — counts  30 

X-axis  calibration  offset  (note  3)  2 

Y-axis  calibration  offset  (note  3)  2 

Z-axis  calibration  offset  (note  3)  2 

calibration  matrix  (3x3)  18 

B-x  deviation,  30  s 

detrend — nT  20 

B-y  devation,  30  s 

detrend — nT  20 

B-x  current  density 

estimate--A/m*m  20 

B-y  curent  density 

estimate--A/m*m  20 

B-x  power  spectal  density 
function  decimated  in  log 

f requency--db  nT  sq/Hz  40 

B-y  power  spectral  density 
function  decimated  in  log 

f requency--db  nT  sq/Hz  40 

Frequencies  of  above  power 

spectral  density  samples--Hz  40 

B  magnitude,  rms  model-data, 

fit--nT  (note  4)  2 

B-sunline  angle,  rms  model-data 
fit--deg  (note  4)  2 

(spares)  12 


i 
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HILAT 


(Cont .  ) 


Number  of 


Particle 

Detector 


Instrument  mode 
(spare) 


Scaled  data,  mode 

dependent--counts  (note  5)  9 

mode  1:  vertical  sensor,  0.25  s  average 
channel  1,  60  samples 
through 

channel  8,  60  samples 


I 

i 


<*3 

ij 

»5 

ij 


mode  2:  vertical  sensor,  0.5  r.  average 
channel  1,  30  samples 
through 

channel  16,  30  samples 

mode  3:  all  sensors,  1.5  s  average 

channel  1,  vertical  sensor,  10  samples 
through 

channel  16,  vertical  sensor,  10  samples 

channel  1,  45  degree  sensor,  10  samples 
through 

channel  16,  45  degree  sensor,  10  samples 

channel  1,  nadir  sensor,  10  samples 
through 

channel  16,  nadir  sensor,  10  samples 

Log  integral  number  flux, energy 
<630  eV,  vertical  sensor — 

e/cm/cm-s-sr  20 

Log  integral  energy  flux,  energy 
<630  eV,  vertical  sensor— 
keV/cm/cm-s-sr  20 

Log  integral  number  fluxenergy 
>1  keV,  vertical  sensor — 

e/cm/cm-s-sr  20 

Log  integral  energy  flux,  energy 
>1  keV,  vertical  sensor— 

keV/cm/cm-s-sr  20 

Integral  energy  flux  energy  <630  eV, 
power  spectral  density  function 
decimated  in  log  frequency 

— db  flux  units  sq/Hz  40 


VI 


HI  LAX 


SCIENCE  SUMMARY  DATA  (Cont . ) 


Parameter 

Number  of 
16-bit  Words 

Particle 

Integral  energy  flux  energy  >1  keV, 

Detector 

power  spectral  density  function 

decimated  in  log  frequency 
— db  flux  units  sq/Hz 

40 

Frequencies  of  above  power 

spectral  density  samples--Hz 

40 

(spares) 

10 

1 DM/HP A 

Ion  density,  IDM--cm“3 

20 

Ion  temperature--degrees  K 

20 

Ion  density,  RPA — cm-3 

20 

(spares) 

20 

Ram  velocity — km/s 

20 

Crosstrack  velocity — km/s 

20 

Vertical  velocity — km/s 

20 

(spares) 

20 

Ion  density  power  spectral  density 
function  decimated  in  log 
f requency--db  e/cm-3  sq/Hz 

16 

Cross-track  velocity  power  spectral 
density  function  decimated  in 
log  frequency — db  km/s  sq/Hz 

40 

Frequencies  of  above  power 

spectral  density  samples--Hz 

40 

(spares) 

40 

Nadir 

3914  A  intensity — rayleighs  20 

Photometers 

6300  A  intensity--rayleighs  20 

Spacecraft 

Rotation  matrix  (3x3),  vehicle  to 

Attitude 

geographic  (note  4) 

18 

B-x  rms  model-data  fit — nT  (note  4) 

2 

B-y  rms  model-data  fit — nT  (note  4) 

2 

B-z  rms  model-data  fit — nT  (note  4) 

2 

hilat 


SCIENCE  SUMMARY  DATA  (Cont.) 


note  1  -  all  parameters  are  HP  32-bit  real  words,  except  for 
magentometer  counts,  which  are  16-bit  integers 


note  2  -  times  are  with  respect  to  the  time/geometry/magnetic  field 
point  and  will  correspond  to  those  in  the  telemetry 
records 

note  3  -  nominal  (February  1984)  values  are  used  if  the  offsets  are 
not  optimized,  e.g.  the  pass  is  sunlit  for  less  than  two 
minutes 

note  4-  dummy  values  (-999)  are  used  if  the  pass  is  sunlit  for  les 
than  two  minutes 


note  5-  the  particle  detector  channel  energies  are  as  follow.'.: 


channel 

1 

.020 

keV 

channel 

2 

.032 

keV 

channel 

3 

.054 

keV 

channel 

4 

.088 

keV 

channel 

5 

.144 

keV 

channel 

6 

.235 

keV 

channel 

7 

.385 

keV 

channel 

8 

.632 

keV 

channel 

9 

.632 

keV 

channel 

10 

1.035 

keV 

channel 

11 

1.700 

keV 

channel 

12 

2.780 

keV 

channel 

13 

4.550 

keV 

channel 

14 

7.450 

keV 

channel 

15 

12.200 

keV 

channel 

16 

20.000 

keV 

APPENDIX  C 


Data  Base  Tapes  9  Track 
6250  bpi 
Labeled 
W-I 

Seneral  Structure: 

Each  tape  has  multiple  passes.  Each  file  has  data  from  one 
station  pass.  Each  file  has  two  preface  records  (preceding  the 
measurement  data) .  There  is  an  EOF  after  each  station  pass  and  a 
iouble  EOF  following  the  data  for  each  experiment. 

Structure  is  as  follows: 

Multiple  passes  of  J  Sensor  (J/SENSOR  Program) 

Double  EOF 

Multiple  passes  of  IDM  density  (TAPE2) 

Double  EOF 

Multiple  passes  of  IDM  drift  velocity  (TAPE3) 

Double  EOF 

Multiple  passes  of  Electron  Sensor  "current"  (TAPE20) 
Double  EOF 

Multiple  passes  of  Electron  Sensor  F1,F2,F3,F4  (TAPE22 ) 
Double  EOF 

Multiple  passes  of  RPA  data  (TAPE21) 

Double  EOF 

Multiple  passes  of  Scintillation  Data  (SCINT.  Program) 
Double  EOF 

Files  may  be  accessed  by  using  a  combination  of  COPYBF  & 

"  PYBR . 

Skip  all  files  for  "Experiments"  using  COPYBF.  Skip  files 
..:*hin  an  experiment  type  by  using  COPYBR. 

For  example,  to  position  tape  at  beginning  of  4th  file  of  RPA 

mat. a  . 

COPYBF, LFN, XX, 5. 

COPYBF, LFN, XX, 3 . 


hilat 

Preface _ Record  1 


HQXd. 

Description 

Format 

1 

Year 

F 

2 

Day  of  year 

F 

3 

Geodetic  latitude  of 

station 

F 

4 

Longitude  of  station 

F 

5-15 

Orbital  elements  used 
computations  for  this 

in  ephemeris 
pass 

F 

EILAT 

Prefaea _ Record _ Z 

The  time  spacing  between  'frames'  is  15  seconds.  There  are 
’ M '  groups  of  45  words  in  this  record  (since  the  maximum  pass 
duration  is  less  than  20  minutes,  this  record  should  never  exceed 
3500  words) .  The  first  word  in  this  record  is  N  (ineger),  the 
number  of  groups  of  45  words.  The  format  of  each  group  in  the 
record  repeated  N-l  times  is  as  follows: 

Word  No.  Description  Format 


1  UT  (seconds)  F 

2  Geodetic  latitude  (deg)  F 

3  Longitude  (+E)  F 

4  Geodetic  altitude  (km)  F 

5  Geocentric  position  X  (km)  F 

6  Geocentric  position  Y  (km)  F 

7  Geocentric  position  Z  (km)  F 

8  Geocentric  velocity  X  (km/s)  F 

9  Geocentric  velocity  Y  (km/s)  F 

10  Geocentric  velocity  Z  (km/s)  F 

11  Orbital  attitude  pitch  F 

12  Orbital  attitude  yaw  F 

13  Orbital  attitude  roll  F 

14  Invariant  magnetic  latitude  F 

15  Magnetic  local  time  (seconds)  F 

16  Magnetic  dip  (degrees)  F 

17  Magnetic  declination  (degrees)  F 

18  Magnetic  strength  -  North  (NT)  F 

19  Magnetic  strength  -  East  (NT)  F 

20  Magnetic  strength  -  Vertical  (NT)  F 

21  Solar  right  ascension  (degrees)  F 

22  Solar  declination  (degrees)  F 

23  Solar  zenith  angle  F 

24  Sun-shade  indicator  F 

25  CGM  latitude  (deg)  F 

26  CGM  longitude  (deg  F 


hilat 

Preface  Record  2  (Cont . ) 


Word  NO.  Description 

27  CGM  local  time  (hrs) 


Format 


28 

MF 

line 

thru 

satellite 

- 

Geod . 

lat . 

(350 

km) 

F 

29 

MF 

line 

thru 

satellite 

- 

Good. 

long 

.  (350  km) 

F 

30 

MF 

line 

thru 

satellite 

- 

Geod . 

alt  . 

(350 

km) 

F 

31 

MF 

line 

thru 

satellite 

- 

Geod . 

SZA 

(350 

km) 

F 

32 

MF  line 
(350  km) 

thru 

satellite 

- 

S/S  indicator 

F 

33 

MF 

line 

thru 

satellite 

- 

Geod . 

lat  . 

(100 

km) 

F 

34 

MF 

line 

thru 

satellite 

- 

Geod . 

long 

.  (100  km) 

F 

35 

MF 

line 

thru 

satellite 

- 

Geod . 

alt  . 

(100 

km) 

F 

36 

MF 

line 

thru 

satellite 

- 

Geod . 

SZA 

(100 

km) 

F 

37 

MF 

line 

thru 

satellite 

- 

Geod . 

S/S 

(100 

km) 

F 

38 

MF  line  thru 
(surface) 

satellite 

- 

Geod . 

lat . 

F 

39 

MF  line  thru 
(surface) 

satellite 

- 

Geod. 

long 

• 

F 

40 

MF 

line 

thru 

satellite 

- 

Geod . 

alt  . 

F 

(surface) 

41  Inv.  Mag.  Lat .  of  F  region  Pen  (350  km) 

42  Magnetic  LT  (hrs)  {0-24} 

43  Vacant 

44  Vacant 

45  Vacant 


F 

p 

F 

F 

F 


HLLAI 

J/SenaQg _ Data  Base 


There  will  be  one  file  per  pass.  Each  file  will  have  a 
header  record  followed  by  a  series  of  data  records. 


Header  Record 

0.1  Word  count  (8)  (I) 

1  Year  F 

2  Day  F 

3  Latitude  of  station  (* (**))  F 

4  Longitude  of  station  (*)  F 

5  Instrument  mode  (at  beginning  of  pass  F 

( 1 . ,  2 .  or  3  . ) 

6  Start  time  of  pass  (UT  seconds)  F 

7  Run  date  of  file  creation  A 

8  Run  time  of  file  creation  A 


(*)  Indicates  information  obtained  from  SRI  tape 

(**)From  AFGL/SUNY  ephemeris 


EILAT 


(Cont  .  ) 


J/Sansor  Data  Base 


Data  Records  (One  second  per  record) 


ft 


K 


,  1 
1 
2 

3 

4 

5 

6 

7 

8 


Word  count  (112) 

Mode  indicator  (1.,  2 
UT  (sec) 

Alt  (km) 

Longitude 
Geodetic  latitude 
Local  time  (sec) 
Geomagnetic  latitude 
Magnetic  LT  (sec) 


or  3 . ) 


(I) 

F 


9 

CGM 

latitude 

10 

CGM 

longitude 

11 

CGM 

LT  (hours  + 

Frac  of 

liDiira.) 

12 

Invariant  latitude 

13 

Mag 

field  Bx  "'l 

) 

14 

Mag 

field  By  ' 

From 

igrf80 

15 

Mag 

field  Bz  J 

i 

16 

Pitch  ''j 

17 

Yaw 

|  From  summary  tape*** 

18 

Roll  J 

19 

x  1 

20 

Y 

21 

Z 

Spacecraft 

position 

22 

X 

f  and  velocity*** 

23 

Y 

24 

Z 

y 

25-35 

Spares 

36 

HV  monitor  (V)  from  most 

recent 

37 

LV  monitor  (V)  from  most 

recent 

From  AFGL/SUNY 
ephemer is 


aiLAT 


Word 

38 


39-86 


87-102 


J/Sensor  Data  Base  /Conti  .  ) 


OfiSGripLion  Format 

Temp  (°C)  (V)  from  most  recent  TLM  F 

readout 

Particle  counts  words  are  the  computed 
counts  from  the  TLM  data  (i.e.  counts 
=  2E (AM+32 ) -33 )  packed  into  15  bit 
words.  Frames  missing  due  to  TLM 
dropout  will  have  1  filled  15  bit 
words  (77777g) .  There  are  192  15  bit 
words  (48  60  bit  words) 

Packed  counts  data  stored  as  follows:  (Z  I 
=  Zenith  detector;  F  =  Forty  deg;  N  = 

Nadir)  -  subscript  indicates  channel 

In  mode  1 : 

Zx,  Z2/ Z3,  . . . , Zg  repeated  24  times 
In  mode  2 : 


Zlr Z2,  .  .  . , Z16  repeated  12  times 

In  mode  3: 


lr Z2>  ■ 

•  •  t  ^16  , 

Zi,  . . 

•  ' z16  > 

ZX,  . 

•  •  t  Z16  / 

ZX,  . 

■  vZ16 

l'F2'  • 

•  •  '  F16  ' 

FX,  .  . 

•  '  F16  ' 

Fl'  • 

•  •  r F16  ’ 

Flr  • 

■  •  '  F1 6 

i,n2,  • 

•  •  >  n16  ' 

NX,  .  . 

•/N16  , 

NX,  . 

•  •  r  n16  > 

NX,  • 

•  -/N16 

(i.e.  4  groups  of  16  zenith  spectra,  followed 
by  4  groups  of  forty  degree  detector  spectra 
followed,  by  4  groups  of  Nadir  detector 
spectra . ) 

The  next  group  of  words  will  be  31  words  I 
per  telemetry  frame  (-1/2  sec)  from 
the  magnetometer  (total  of  62  words) 
packed  into  15  bit  words.  (The  30  LSBS 
of  word  102  will  be  vacant) 


103 


ID  word  from  even  block 


I 


EILAT 


J/Sensor  Data  Base  (Cont . ) 


Description  Format 

Time  of  telemetry  frame  for  which  Line  F 

No.  =  6. 

Packed  photometer  data  (56  bits  of  I 

information  -right  adjusted  -  the  6 
photometer  words  plus  housekeeping  word 
41  which  contains  photometer  on/off 
discretes . ) 


65 


yMV'VS1 


1 


Year  (last  2  digits  of  19xx) 
Day  of  year 
UT  (seconds) 


2 

3 

4  IDM  frame  counter 

5  Ion  density  (cm3) 

Dummy  fill  value  is  -999999. 

IDM  DATA  BASK  l  Drift 

(20  words/record) 

Word  Description 

1  Year  (last  2  digits  of  19xx) 

2  Day  (day  of  year) 

3  UT  (seconds) 

4  IDM  frame  counter 

3-20  16  velocity  values  [order  is 

frame  counter] 


Dummv  fill  value  is  -999999. 


HI  LAX 

ELECTRON  -SENSOR  DATA  BASE  (Cu rrftnM 


(20  words/record) 
Mold  Descript  j on 


1 

2 

3 

4 

5-20 


Year  (last  2  digits  of  19xx) 

Day  (day  of  year) 

UT  (seconds) 

Electron  sensor  frame  counter 
16  values  of  electron  sensor  tlm.  volts 


HOXiHsiL 


F 

T 

P 


£ 


Dummy  fill  value  is  -999999. 

ELECTRON  SENSOR  DATA  BASE  fPj^F^ ,  F2-Fjl_ 

(8  words/record) 

Word  Description 


1 

2 

3 

4 

5 

6 

7 

8 


Year  (last  2  digits  of  19xx 
Day  (day  of  year) 

UT  (seconds) 

Electron  sensor  frame  counter 
Fx  (volts) 

F2  (Volts) 

F3  (volts) 

F4  (volts) 


Format 


I 

I 

F 

I 

F 

p 

p 

F 


F.mmy  fill  value  is  -999999, 


67 


hilat 


(35  words/record) 
Word  Descript 


19-35 


Year  (last  2  digits  of  19xx)  I 

Day  (day  of  year)  I 

UT  (seconds)  F 

RP A  frame  counter  I 

Velocity  (km/s)  F 

Sensor  potential  F 

Vehicle  potential  F 

Satlelite  potential  w.r.t.  instrument  bias  F 

lion  temperature  (deg  K)  F 

Ion  density  1  (cm3)  F 

Mass  1  (kg)  F 

Ion  density  2  (cm3)  F 

Mans  2  (kg)  F 

Drift  velocity  (m/s)  F 

Mass  in  AMU  of  light  component  (H+=l,  HE+=4)  I 

Total  ion  density  (cm3)  F 

Est .  variance  of  the  linear  ion  current  F 

Est .  variance  of  the  logarithmic  ion  current  F 

FA ( I , K)  = Array  containing  final  results  F 

(1=1, 11  and  K=l, 3) 

FA ( I , 1 )  Contains  values  of  the  "Best  Solution" 

with  respect  to  the  minimum  of  the  linear 
variance,  where  only  one  parameter-position 
of  the  second  characteristic  slope  is 
allowed  to  vary. 

FA ( I , K)  Where  K=2,3,  contain  values  of  the  limiting 

solutions  to  the  left/right  of  the  "Best" 

position  of  the  second  characteristic  slope 

FA ( 1 , K)  Est.  var.  of  the  lin.  ion  current 

FA(2,K)  Est.  var.  of  the  log.  ion  current 

FA(3,K)  VHALF ( 2 )  position  of  the  second  char,  slope 


FA (4 , K)  Slope  at  VHALF ( 2 ) 


m 


HI  LAX 

RPA  DATA  BASE  (Cont  .  ) 


(35  words/record) 

Wora  Description  Format 

19-35  FA ( 5, K)  Ion  current  at  VHALF ( 2 ) 

FA ( 6, K)  Ion  density  1  (cm3) 

FA  ( 7 ,  K)  Ion  density  2  (cm3) 

FA(8,K)  Ion  temp,  (deg  K) 

FA(9,K)  Drift  velocity  normal  to  RPA  apert .  (m/s) 

FA ( 10, K)  Satellite  potential  w.r.t.  instrument  bias 
(volts) 

FA ( 1 1 , K)  Mode,  1=1  ion  case,  2=2  ion  case 
For  no  solution 

Dummy  fill  value  is  9.9  x  109  or  9.9  x  103 


HI  LAI 


Desc . 

Year  (I)  e.g.  83,84,85, 

Day  of  year  (I) 

Station  code  (I)  (l=Sondre;  2=Tromso,  3=Churchill; 
4=Rover) 

UT  (seconds) 

Words  5-69  from  TIME/GEOM/MF  Record 

Desc. 


-  station  Geodetic  latitude  -  degree  +N 

Longitude  -  degree  +E 

Geodetic  altitude  -  km 

Azimuth  angle  to  satellite  -  degree 

Elevation  angle  to  satellite 
-  degree 


-  satellite 


Geodetic  latitude  -  degree  +N 

Longitude  -  degree  +E 

Geodetic  altitude  -  km 

Orbital  attitude  pitch  -  degree 

Orbital  attitude  yaw  -  degree 

Orbital  attitude  roll  -  degree 

Invariant  magnetic  latitude, 

APL  eccentric  dipole  -  degree 

Magnetic  local  time, 

APL  eccentric  dipole  -  0  seconds 

Magnetic  dip,  IGRF80  -  degree 
Magnetic  declination,  IGRF80 
-  degree 

Solar  right  ascension  -  degree 
Solar  declination  -  degree 
Solar  zenith  angle  -  degree 


'V,v!S' - 


hilat 


-  magnetic 
field  line 

through 

satellite 

(IGRF80) 


Geodetic  latitude  (350  km) 

-  degree  +N 

Longitude  (350  km)  -  degree  +E 
Geodetic  altitude  (350  km)  -  km 

Geodetic  latitude  (100  km) 

-  degree  +N 

Longitude  (100  km)  -  degree  +E 
Geodetic  altitude  (100  km)  -  km 


-  F-region  Geodetic  latitude  -  degree  -t-N 
penetration  Longitude  -  degree  +E 
(350  km)  Geodetic  altitude  --  km 

Invariant  magnetic  latitude, 

APL  eccentric  dipole  -  degree 

Magnetic  local  time, 

APL  eccentric  dipole  -  seconds 

Magnetic  dip,  IGRF80  -  degree 

Magnetic  decliaton,  IGRF80  -  degree 

Magnetic  strength  north,  IGRF80  -  n 

Magnetic  strength  east,  IGRF80  -  mg 

Magnetic  strength  vertical,  IGRF8q 

-  mg 

Propagation  zenith  angle  --  degree 

Propagation  magnetic  azimuth 

-  degree 

Off-magnetic  field  propagation 
angle,  IGRF80  -  degree 

Off-magnetic  meridian  propagation 
angle,  IGRF80  -  degree 

Off-magnetic  L-shell  propagation 
angle,  IGRF80  -  degree 

Reduced  propagation  range  -  km 

Magnetic  penetration 
velocity  X  -  km/s 

Magnetic  penetration 
velocity  Y  -  k/S 


m 


ILLATION  STATISTICS  DATA  BASE  (Cont  .  ) 
Words  5-69  from  TIME/GEOM/MF  Record 


isc . 

Magnetic  penetration 
velocity  Z  -  km/ s 

E-region  Geodetic  latitude  -  degree  +N 
penetration  Longitude  -  degree  +  E 
(100  km)  Geodetic  altitude  -  km 

Invariant  magnetic  latitude, 

APL  eccentric  dipole  -  degree 

Magnetic  local  time, 

APL  eccentric  dipole  -  seconds 

Magnetic  dip,  IGRF80  -  degree 

Magnetic  declination,  IGRF80 

-  degree 

Magnetic  strength  north,  IGRF80  - 

Magnetic  strength  east,  IGRF80  -  i 

Magnetic  strength  vertical,  IGRF8 

-  mg 

Propagation  zenith  angle  -  degree 

Propagation  magnetic  azimuth 

-  degree 

Off-magnetic  field  propagation 
angle,  IGRF80  -  degree 

Off-magnetic  meridian  propagation 
angle,  IGRF80  -  degree 

Off-magnetic  L-shell  propagation 
angle,  IGRF80  -  degree 

Reduced  propagation  range  -  km 

Magnetic  penetration 
velocity  X  -  km/s 

Magnetic  penetration 
velocity  Y  -  km/s 

Magnetic  penetration 
velocity  Z  -  km/s 


hilat 


RW.wf-V  *  ■■  1 1  n  I  \  i 


80 

81 

82-84 


(Cont . ) 

Words  5-69  from  TIME/GEOM/MF  Record 


>rd  No. 

Desc . 

67 

-  Magnetic 

Geodetic  latitude 

(100  km) 

field  line 

-  degree  +N 

through 

68 

F-region 

Longitude  (100  km) 

-  degree  +E 

penetration 

69 

Geodetic  altitude 

(100  km)  -  km 

Words  70-84  from  Scintillation  Records 
Desc . 


-137  MHz 


413  MHz 
antenna  1 


-  1239  MHz 


Mean  signal  level  -  db 
Intensity  scintillation  index  S4 
Intensity  fade  period  -  second 

Intensity  decorrelation  time 

-  second 

Standard  deviation  of  phase  -  rad 

Mean  signal  level  -  db 
Intensity  scintillation  index  S4 
Intensity  fade  period  -  second 

Intensity  decorrelation  time 

-  second 

Standard  deviation  of  phase  -  rad 

Mean  signal  level  -  db 

Intensity  scintillation  index  S4 

Total  electron  content  along 
propagation  ray  path,  three 
equispaced  samples  -- 
e  /m  sq 


EILAT 

SCINTILLATION  SSTATISTICS  DATA  BASE  (Cont.) 
Words  85-204  from  SCI/SUM  Records 


Word  No.  Desc 


85-94 

95-104 

105-114 

: 15-124 

125-134 


-  Particle 
Detector 


-  IDM/RPA 


Log  integral  number  flux, 
energy  <630  ev,  vertical  sensor  — 
e/cm/cm-s-sr,  ten  samples 

Log  integral  energy  flux, 
energy  <630  ev,  vertical  sensor  -- 
kev/cm/cm-s-sr,  ten  samples 

Log  integral  number  flux, 
energy  >1  kev,  vertical  sensor  -- 
e/cm/cm-s-sr,  ten  samples 

Log  integral  energy  flux, 
energy  >1  kev,  vertical  sensor  -- 
kev/cm/cm-s-sr,  ten  samples 

Ion  density,  IDM  --  cm-3,  ten 

samples 


135-144 

115-154 


Ram  velocity  —  km/s,  ten  samples 

Crosstrack  velocity  --  km/s,  ten 
samples 


155-164 


Vertical  velocity  —  km/s,  ten 
samples 


165-184 

185-194 

195-204 


Frequencies  of  cross  track  velocity 
power  spectral  density  samples  -- 
Hz,  twenty  samples 

Nadir  3914  A  intensity-rayleights,  10 

photometers  samples 

6300  A  intensity  -  rayleighs,  10 
samples 


As  done  for  previous  'scintillation  data  base'* 
Word  No.  Desc . 


205  Slope  (p) 

206  RMS 


^intensity 


137  MHz 


207 


TI  (1  Hz  value) 


74 


HILAT 


(Cont 


phase 


137  MHz 


Hasc . 

Slope  1  (PI) 
RMS  1 

Slope  2  (P2 ) 


211  RMS  2 

212  TP  (1  Hz  value)! 

213  Slope  (P)  \  A 

214  RMS  L intensity 

215  TI  (1  Hz  value) J 

216  Slope  1  (PI)  413  MHz 

217  RMS  1  ]  f 

218  Slope  2  (P2)  Lphase 

219  RMS  2  1 

220  TP  (1  Hz  value)) 

Fits  and  integrations  to  science  summary  data  parameters 


221  Slope/10  (P) 

222  1  Hz  value  (T)  from  By  fit  1 

22  3  RMS  value  (RMS) 

224  Area  (A) 

225-228  As  above  for  By  fit  2 

229-232  As  above  for  integral  energy  flux  <630  eV  fit: 

233-236  As  above  for  integral  energy  flux  <630  eV,  fit  I 

237-240  As  above  for  integral  energy  flux  >1  kev,  fit  1 

241-244  As  above  for  integral  energy  flux  >1  kev,  fit  7 

245-248  As  above  for  cross  track  velocity,  fit  1 

249-252  As  above  for  cross  track  velocity,  fit  1 

253-256  As  above  for  ion  density,  1  fit 

257  x 

258  y  satellite  velocity  (km/sec) 


262-270 


Range  to  satellite  (km) 

J  sensor  instrument  mode 
Vacant 


M'' \v£V\ r >■& 


HILAT 


(24  Parameters) 

Dimension  STMEAN (12, 12) , ST10P (12, 12) , ST50P (12, 12) , 

ST90P (12, 12) , STCNTS (12, 12) 

DO  10  KP=1 , 2 

DO  1  IREAD=1, 24 

1  READ (1) STMEAN, ST10P, ST50P , ST90P , STCNTS 

READ (1)  NFILES , NR1 , NR2 , NR3 , NR4 , NR5, NR 6, NR7 , NR8, ISTNCD, 
IDBEG, IDEND, NR 9, NR 10 

10  Continue 

Five  arrays  are  read  for  each  of  24  parameters: 

STMEAN  is  mean  value 

ST10P  is  10th  percentile  value 

ST50P  is  50th  percentile  value 

ST90P  is  90th  percentile  value 

STCNTS  is  number  of  occurrences  in  each  bin 


Lastly,  a  preface  record  where  NFILES  = 
files  (i.e.  revs)  in  data  set. 

total 

number 

of 

NR1 

-  Input  data  recs  used  after  elevation, 
exclusion  for  parameters  19  &  20 

l'nv . 

Let .  & 

KP 

NR8 

-  Input  data  recs  used  after  elevation. 

Inv . 

Let .  & 

Kp 

exclusion  for  parameters  1  to  18 

NR2  -  Rees  used  after  exclusion  for  scintillation  freq's 
all  zero 

NR3  -  Rees  remaining  after  S4  exclusion 

NR 4  -  Rees  remaining  after  exclusion 

NR5  -  Rees  remaining  after  10  ion  density  values  outside 
range  103  to  106 

NR6  -  Rees  remaining  after  avg .  energy  exclusion  for  Mode 
*2  or  3 

NR7  -  Rees  excluded  for  avg.  energy  denominator  <103 

ISTNCD  -  Station  code  (l=Sondre;  2=Tromso;  3=Churchill) 

TDBEG  -  Start  day  of  data  set  (e.g.  83317) 


IDEND  -  End  day  of  data  set  (e.g.  84031) 

NR9  -  Rees  excluded  from  low  energy  of  flux  for 
denominator  <1000 

NR10  -  Rees  excluded  from  high  energy  &  flux  for 
denominator  <1000 

This  sequence  is  done  twice  for  2  Kp  values:  Kp  <3.5  and  Kp 
>3.5. 


For  all  arrays,  the  values  are  stored  as  12x12,  with  rows 
1->12  being  Magnetic  Local  Time,  starting  at  23  hours, 
incrementing  in  2  hour  bins.  Columns  1-»12  are  the  Invariant 
Latitude  bins,  starting  at  50°,  incrementing  in  2.5°  bins. 

In  the  case  of  1  hour  magnetic  local  time  bins,  for  all 
arrays,  the  values  are  stored  as  24x12,  with  rows  1-24  being 
magnetic  local  time,  starting  at  22.5  hours,  incrementing  in 
1  hour  bins.  Columns  1-12  are  the  invariant  latitude  bins, 
starting  at  50°,  incrementing  in  2.5°  bins. 


HILAT  SCINTILLATION  DATA  BASE  STATISTICS  FILE  FORMAT  (Cont . ) 

(24  Parameters) 


The  names  of  the  24  parameters  are: 
No .  Name 


1 

2 
3 

A 

*1 

5 

6 

7 

8 
9 

10 
1 1 
12 
:  3 
14 
1  5 

:  6 
17 
19 

1  9 

r\ 

i.  i 

2  2 

; 

2  4 


Mean  signal,  137  MHz 

Scintillation  index  intensity,  137  MHz 
Fade  period  intensity,  137  MHz 
Decorrelation  time  intensity,  137  MHz 
Phase  std.  deviation,  137  MHz 
Mean  signal,  413  MHz 

Scintillation  index  intensity,  413  MHz 
Fade  period  intensity,  413  MHz 
Decorrelation  time  intensity,  413  MHz 
Phase  std.  deviation,  413  MHz 
Mean  signal  1239  MHz 
TEC 

Phase  slope  1,  137  MHz 

Phase  slope  2,  137  MHz 
Phase  Tp,  137  MHz 

Phase  slope  1,  413  MHz 

Phase  slope  2,  413  MHz 
Phase  Tp,  413  MHz 

Average  ion  density 
Average  energy-particle  detector 
Average  energy  -  particle  detector  (low) 
Energy  flux  -  particle  detector  (low) 
Average  energy  -  particle  detector  (high) 
Energy  flux  -  particle  detector  (high) 
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.VJmLVjVJA’a'j  >jA 


Dimension  STMEAN ( 12 , 12 )  ,  ST10P  ( 12,  12) ,  ST50P ( 12 , 12 ) , 
ST90P  (12, 12) , 

STCNTS  (12, 12) 

DO  10  KP=1 , 2 

DO  1  IREAD=1 , 2  7 

1  READ ( 1 )  STMEAN, ST10P, ST50P,  ST90P,  STCNTS 

READ (1)  NFILES, (NR(I) , 1=1, 10) , ISTNCD, IDDBEG, IDEND 


10  Continue 


Five  arrays  are  read  for  each  of  27  parameters: 


STMEAN 

ST10P 

ST50P 

ST90P 

STCNTS 


is  mean  value 

is  10th  percentile  value 

is  50th  percentile  value 

is  90th  percentile  value 

is  number  of  occurrences  in  each  bin 


Lastly,  a  preface  record  where: 

NFILES  =  Total  number  of  files  (i.e.  revs) 

in  data  set 


(NR ( I )  ,  1  =  1 ,  10 )  =  1  total  number  of  input  data 

base  records  read. 

2-10)  records  used  for  each  set  ot 
3  parameters,  respectively. 

ISTNCD  =  Station  code  (l=Sondre;  2=Tromso; 

3=Churchill) 

=  Start  day  of  data  set  (e.g.  84032) 
=  End  day  of  data  set  (e.g.  84121) 
done  twice  for  2  Kp  values:  Kp  <3.5  and  K . 


IDBEG 

IDEND 


it  rays,  the  values  are  stored  as  12x12,  with  rows  1- 
-  ;  -asr.^-ic  local  time,  starting  at  23  hours, 

.'  i  in  2  hour  bins.  Columns  1-12  are  the  invariant 
•  .  n  .  starting  at  50°,  incrementing  in  2.5°  bins. 

>  ■■  f  1  hour  magnetic  local  time  bins,  for  all 
,  •  •  •*  vi  lues  are  stored  as  24x12,  with  rows  1-24  be'ri 
•  :me,  starting  at  22.5  hours,  incremen- no  ,  : 
Columns  1-12  are  the  invariant  latitude  hi:.:, 

°,  incrementing  in  2.5°  bins. 


1 


HILAT  SCIENCE  SUMMARY  DATA  BASE  STATISTICS  FTT.F.  format 

(Cont  . ) 

(27  Parameters) 


The  names  of  the  27  parameters  a 


re  : 


Name 

i 

By  slope  1 

2 

By  T1 

3 

By  area  1 

4 

By  slope  2 

5 

By  T2 

6 

By  area  2 

7 

Loen  slope  1 

9 

Loen  T1 

9 

Loen  area  1 

10 

Loen  slope  2 

1 1 

Loen  T2 

12 

Loen  area  2 

13 

Hien  slope  1 

14 

Hien  T1 

15 

Hien  area  1 

16 

Hien  slope  2 

17 

Hien  T2 

1  8 

Hien  area  2 

19 

CRVEL  slope  1 

CRVEL  T 1 

■’ 

CRVEL  area  1 

.  .. 

CRVEL  slope  2 

2  3 

CRVEL  T2 

/  4 

CRVEL  area  2 

.  - 

Ion  den  slope  1 

Ion  den  T1 

<  7 

Ion  den  area  10 
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hilat 

PHOTOMETER  DATA  BASE 


Record  1 
Record  2 
Record  3 

MQXd— Mn,. 

1-4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 
23 
29 

30-4  0 
4  1-87 
88-95 
95 

96-100 


Input  pref  record  1  from  input  (summary  tape) 
Input  ephemeris  rec  from  input  (summary  tape) 
100  words  as  follows: 

Not  used 
UT  (secs) 

Geod.  altitude 
Geod.  longitude 
not  used 
Geod.  latitude 
not  used 
CGM  latitude 
CGM  longitude 
CGM  local  time 
Solar  zenith  angle 
Sun-shade  indicator 
Geocentric  position  -  X 
Geocentric  position  -  Y 
Geocentric  position  -  Z 
Geocentric  velocity  -  X 
Geocentric  velocity  -  Y 
Geocentric  velocity  -  Z 
Orbital  attitude  pitch 
Orbital  attitude  yaw 
Orbital  attitude  roll 
Inv.  Mag.  latitude 
Mag.  local  time 
Mag.  strength  -  North 
Mag.  strength  -  East 
Mag.  strength  -  Vertical 
Not  used 

Magnetometer  (6  frames) 

Photometer  &  Housekeeping  words 
ID 

Not  used 
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HILAT 


AIM  Data  Base 

There  will  be  one  file  per  station  pass.  Each  file  will 
have  a  2  header  records  followed  by  a  series  of  data  records . 


Header  Record  1 

1  Year  F 

2  Day  F 

3  Latitude  of  station  F 

4  Longitude  of  station  F 

5  AIM  mode  at  beginning  of  pass  F 

1.  =  Imaging; 

2.  =  Spectrometer; 

3.  =  photometer 


6  Start  time  of  pass  -  UT  seconds  F 

7  Run  date  of  file  creation  A 

8  Run  time  of  file  creation  A 

9  Wavelength  (if  in  imaging  or  photometer  F 

mode  -  if  spect  mode  value  =  -1.) 


Header  Record  2 

Epnemeris  data  REC  (45xN) 

Same  as  Preface  REC  2  under  Input,  part  VIII 


Data  Records 

Each  data  record  is  constructed  from  6  telemetry  frames 
(~3  seconds) .  The  record  is  structured  such  that  the  first 
AIM  word  in  a  data  record  is  the  first  word  of  a  mode  cycle. 
Data  from  the  imager  is  packed  into  15  bit  words.  Telemetry 
words  missing  due  to  dropout  will  be  1  filled  (i.e.  77777g) . 


ft 


m 


fij 


2 

3 


5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


Spect  rometer ; 
2  .  = 


Mode  indicator  (1.  =  Imaging;  2. 

3.  =  photometer) 

Code  (0.  =  Normal;  1.  =  Electrical  Test; 

Optical  Test) 

Status  (0.  =  Normal;  1.  =  Same  words  filled  due  to 
TLM  dropout) 

Wavelength  (meaningful  only  in  imaging  & 
photometer  mode 

UT  (sec) 

Altitude  (km) 

Longitude  (+E) 

blank 

Geodetic  latitude 
blank 

Corrected  Geomagnetic  latitude 

Corrected  Geomagnetic  longitude 

Corrected  Geomagnetic  local  time 
(sec) 

Solar  zenith  angle 

Sun/shade  indicator  (-1.  =  shade; 

+1.  =  sun) 

X  ^ 

Y  f  position  vector 
Z 
X 

Y 
Z 

Pitch 
Yaw 
Roll 

Invariant  mag.  lat 
Magnetic  LT 
bN 

bE 
Bv 


from  Preface 

► 

Rec  2-Ephemeri: 


VI 


VtV 

I 

Xt* 

;v;t 


M 

$ 

s?«*% 


$ 

v® 

wi 

>$ 


$ 

Sf. 

tV;4 


33-40 

41-87 


88-171 

172-180 


181-190 


Sun  sensor  indicated 


0  =  all  s/s  off; 

1  =  photo  A  on; 

2  =  photo  B  on 

3  =  A  &  B  on ; 

4  =  AIM  on; 

5  =  A  &  AIM  on; 

6  =  B  &  AIM  on; 

7  -  all  on) 


Line  counter 


Prescalar  flag  (0.  =  no  prescalar;  1.  = 
prescalar) 

Vacant 

Magnetometer  data  (31  words  per  telemetry 
frame  x  6  frames  are  packed  into  15  bit  words. 

This 

requires  46.5  60-bit  words;  the  30  LSBS  of  word  47 


vacant  magnetometer  data 
Pixel  data  (336  15  bit  words 


84  CDC  words) 


The  remaining  words  (15  bits/word)  in  the  record 
are  in  the  sequence  below  and  result  in  9  CDC 
words : 

Fixed  photometer  data  (6  15-bit  wds) 

X  position 
(l)Nadir  point  (1) 

Spectrometer  Status  words  (6) 

AIM  housekeeping  (16) 

ID  words  from  scan  (6) 

Vacant  (spares) 


APPENDIX  D 

AIRS  DATA  BASE  FORMAT 


HXLAX 

AIRS  Data  Base  Format 

The  AIRS  data  base  is  created  with  one  station  pass  per 
file.  The  data  base  is  archived  on  9-Track  digital  tape  at  a 
density  of  6250cpi.  Each  tape  contains  the  data  from  at 
least  one  Summary  Tape. 

Data  base  files  will  consist  of  2  preface  (or  header) 
records  followed  by  a  series  of  AIRS  data  base  records. 

PREFACE  RECORD  1 

This  record  consists  of  15  words. 


Word  No.  Description 


1  Year  (F) 

2  Day  of  year  (F) 

3  Latitude  of  station  (F) 

4  Longitude  of  station  (F) 

5  AIRS  operating  mode  at  beginning  of  the 
pass  (F) 

l.=Imaging,  2 . =Spectrometer ,  3 . =Photometer 

6  Start  time  of  pass  -UT  seconds  (F) 

7  Run  date  of  file  creation  (A) 

8  Run  time  of  file  creation  (A) 


9 

Wavelength 

mode) 

of 

Det 

1 

(-1  • 

if  in  spectrometer 

10 

Wavelength 

of 

Det 

2 

( 

tl  ll  ll 

1 1 

Wavelength 

of 

Det 

3 

( 

it  it  ll 

12 

3-15 

Wave  length 

Vacant 

of 

Det 

4 

( 

ll  ll  n 

;  :■  Er  ACE  RECORD  2 

This  record  contains  satellite  ephemeris,  magnetic  and 
<»*.t itude  data  extracted  from  the  Summary  Tape.  Data  stored  in 
• h i s  record  is  written  in  block  form  to  allow  for 
interpolation  of  parameters  not  stored  on  data  base  records. 
The  time  spacing  between  blocks  is  15  seconds  (as  is  on  the 
nummary  Tape)  .  There  are  45  parameters  contained  in  each 
j  iOck.  A  leading  integer  word  is  provided  to  indicate  the 
number  of  blocks  in  the  record  (since  record  size  will  vary 
impending  on  the  time  duration  of  the  pass) .  No  pass  should 
•■xceed  30  minutes  and  thus  the  record  will  not  be  larger  than 
5401  words.  All  words  are  60  bit  floating  point. 
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RLLAI 

AIRS  Data  Base  Format  (Cont.) 

A  READ  statement  of  the  form 

READ (x)  N,  ((D(I,J),  1=1,45),  J=1,N) 
should  be  used  to  input  the  entire  record. 


The  word  order  of  the  45  words  in  each  block  is  as 
follows : 


Word  No.  Description 


1 

2 

3 

4 

5 


UT  -  seconds 
Geodetic  latitude  -  deg 
Longitude  (+E)  -  deg 
Geodetic  altitude  -  km 
Geocentric  position  X  -  km 


r 

II 

II 

Y  - 

km 

7 

II 

II 

Z  - 

km 

8 

Geocentric  velocity  X  - 

km/ sec 

9 

II 

II 

Y  - 

km/sec 

10 

It 

II 

Z  - 

km/sec 

11 

Pitch  -  deg 

12 

Yaw  -  deg 

13 

Roll  -  deg 

14 

Invariant 
dipole) - 

mag . 
deg 

latitude 

(APL  eccentric 

15 

Magnetic 

-sec 

local 

time  (APL 

eccentric  dipoL 

16 

Magnetic 

dip  - 

deg  (IGRF  85) 

17 

Magnetic 

declination  -  deg  (IGRF  85) 

18 

Magnetic 

strength  North 

-  mG  (IGRF  85) 

19 

II 

II 

East  - 

mG  (IGRF  85) 

20 

It 

ll 

Vertical  -  mG  (IGPF  8 

21 

Solar  right  ascension  - 

deg 

22 

Solar  declination  -  deg 

23 

Solar  zenith  angle  -  deg 
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hilat 


AIRS  Data  Base  Format  (Cont 


Word  No. 

Description 

! 

24 

Sun/shade  angle  -  deg 

25 

Invariant  mag.  lat 

(Gustaf sson 

model!  - 

deg 

M 

26 

Magnetic  longitude 

(Gustaf f son 

model)  - 

deg 

PS 

27 

Magnetic  local  time 

(Gustaf sson 

.  model) 

-hrs 

Words  28  through  40  contain  information  on  the  magnetic  field 
line  through  the  satellite  projected  to  350  km,  100km,  and 
the  earth's  surface. 


350km  geodetic  latitude  -  deg 

"  longitude  -  deg 

"  altitude  -  km 

"  solar  zenith  angle 

"  Sun/shade  angle  -  deg 

100km  geodetic  latitude  -  deg 

"  longitude  -  deg 

"  altitude  -  km 

"  solar  zenith  angle  -  deg 

"  Sun/shade  angle  -  deg 

Surface  geodetic  latitude  -  deg 
"  longitude  -  deg 

"  altitude  -  km 

Inv.  Mag  Lat  of  F-region  penetration  -  deg 

Magnetic  local  time  of  F-region  penetration 
-  hr s 

Code  word  (0.  =  no  science  summary  record  on 
Summary  Tape;  1.  =  science  summary  record  on 
tape) 

Vacant 

Vacant 
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HIL&X 


AIRS  Data  Base  Format  (Cont.) 


AIRS  DATA  BASE  RECORDS 

There  are  419  words  per  record. 

Word  No.  Description 


1 

3 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


Mode  indicator  { 1 . =Imaging;  2 . =Spectrometer ; 
3 . =Photometer ) 


Test  status  (l.=Normal;  2.=Dark  shutter; 
3.=Optical  test;  4.=Extended  dark  shutter) 


TLM  status  (O.=normal;  l.=some  words  l’s  filled) 


Detector  1  wavelength  (-1.  if  spectrometer  mode) 
Detector  2  wavelength  (  "  "  "  "  ) 

Detector  3  wavelength  (  "  "  "  "  ) 

Detector  4  wavelength  (  "  "  "  "  ) 

UT  -  sec 
Altitude  -  km 
Longitude  -  deg 
Vacant 

Geodetic  latitude  -  deg 
Vacant 

Invariant  mag.  lat  (Gustafsson  model)  -  deg 
Magnetic  longitude  (Gustafsson  model)  -  deg 
Magnetic  local  time  (Gustafsson  model)  -hrs 
Solar  zenith  angle 
Sun/shade  angle  -  deg 
X  position  -  km 

Y  position  -  km 
Z  position  -  km 

X  velocity  -  km/sec 

Y  velocity  -  km/sec 
Z  velocity  -  km/sec 
Pitch  -  deg 

Yaw  -  deg 
Roll  -  deg 


ck] 


HILAT 

AIRS  Data  Base  Format  (Cont.) 


Word  No. 

28 

deg 

29 

30 

31 

32 

33 

34 

35-40 

41-87 


88-169 

170-251 

252-333 

334-415 


416-419 


Description 

Invariant  mag.  latitude  (APL  eccentric  dipole) - 


Magnetic  local  time  (APL  eccentric  dipole) -sec 

B-North  -  mG 
B-East  -  mG 
B-Vertical  -  mG 
SDF  frame  count 
AIRS  line  counter 
Vacant 

Magnetometer  data  stored  in  15  bit  words. 

(There  are  31  15  bit  magnetometer  words  for  each 
SDF  frame.  There  are  10  words  for  the  X,  Y  and  Z 
magnetometers  and  one  flag  word.  The  groups  of  10 
words  represent  one  B-field  measurement  and  9 
delta  B  values.  Six  such  frames  (accumulated 
over  the  AIRS  scan)  will  be  stored.  The  30  LSB ' s 
of  word  87  will  be  vacant. 

Detector  1  pixel  data 

Detector  2  pixel  data 

Detector  3  pixel  data 

Detector  4  pixel  data 

(The  telemetry  counts  for  the  326  pixel  words  from 
each  detector  will  be  decompressed  and  stored  in 
15  bit  words.  This  results  in  82  CDC  words  per 
detector.  The  30  LBS's  of  the  82nd  word  for  each 
detector  will  be  vacant.) 

AIRS  status  and  housekeeping  words. 

(The  16  telemetry  words  (206-221)  from  line  count 
6  will  be  retained  in  counts  and  stored  into  15 
bit  words. 
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APPENDIX  E 


V* 


ft 182  999  ANALVSIS  SVSTEHS  FOR  RIB  FORCE  HISSIONSCU)  BOSTON  COLL 

CHESTNUT  HILL  HA  SPACE  DATA  ANAL VS IS  LAB 
0  E  OELOREV  ET  AL  28  FEB  87  BC-SDAL-87-1 
UNCLASSIFIED  AFGL-TR-87-8188  F19628-8T-C-8188 


F/G  14/2 


N)  f-J 


(Cont . ) 


NP5  1982  DENSITY  DATA  BASE  TAPES  (20/15  FILTER) 
NINE  TRACK  6250  BP  I  UNLABELLED  (W-I) 

(USE  SKIPF  OR  COPYBF) 


FILE 

CC3257 

CC1158 

CC328 

1 

- 

141 

149 

-> 

134 

142 

151 

3 

135 

143 

152 

4 

136 

144 

153 

5 

137 

145 

154 

6 

138 

147 

157 

7 

139 

148 

158 

8 

140 

155 

159 

9 

146 

156 

160 

10 

162 

:  l 

163 

12 

164 

13 

165 

14 

174 

15 

175 

16 

176 

17 

18 
19 

0 

1 


t  ILE 

CC4072 

CC4472 

CC4473 

1 

235 

259 

271 

-> 

L. 

240 

260 

272 

3 

241 

261 

273 

4 

244 

262 

274 

5 

245 

263 

275 

6 

246 

264 

276 

7 

247 

265 

277 

8 

252 

266 

278 

9 

253 

267 

279 

10 
1 1 
12 


CC3322 

CC1397 

CC4071 

CC1508 

161 

185 

217 

220 

166 

186 

226 

227 

167 

187 

228 

235 

168 

188 

229 

236 

169 

189 

230 

237 

180 

191 

231 

248 

184 

192 

232 

249 

190 

193 

233 

250 

199 

194 

234 

251 

200 

195 

258 

201 

196 

268 

202 

^  97 

269 

214 

i.  98 

270 

215 

206 

283 

216 

207 

284 

208 

285 

218 

290 

219 

291 

305 

306 

307 

CC4498 

CC4501 

CC4503 

280 

301 

313 

281 

302 

321 

282 

303 

322 

287 

304 

323 

288 

308 

324 

295 

309 

325 

296 

310 

326 

- 

311 

327 

312 

328 

329 

330 

331 
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ACCELEROMETER 

DATA  BASE  TAPES  AND  FORMATS  (Cont . ) 

NP 6  1983-1984  DENSITY  DATA  BASE  TAPES 
NINE  TRACK  6250  BPI  UNLABELLED  (W-I) 


(USE 

SKIPF  OR 

COPYBF) 

FILE 

CC4221 

CC3206 

CC4294 

CC3209 

CC3228 

CC3232 

CC4327 

CC4353 

1 

201 

214 

265 

274 

332 

1 

37 

58 

2 

202 

215 

266 

275 

333 

2 

38 

59 

3 

203 

216 

273 

282 

334 

3 

39 

60 

4 

204 

228 

276 

283 

338 

4 

40 

61 

5 

205 

229 

277 

284 

339 

5 

42 

62 

6 

206 

230 

278 

294 

340 

6 

43 

63 

7 

207 

234 

285 

300 

341 

7 

44 

64 

8 

213 

235 

289 

304 

342 

8 

45 

65 

9 

217 

236 

290 

309 

343 

18 

46 

66 

10 

218 

237 

291 

315 

344 

- 

47 

67 

11 

219 

238 

292 

316 

345 

- 

48 

68 

12 

220 

239 

293 

317 

346 

29 

49 

69 

13 

221 

248 

301 

318 

347 

30 

50 

70 

14 

222 

249 

302 

319 

354 

31 

51 

71 

15 

247 

250 

303 

320 

355 

32 

52 

72 

16 

257 

258 

310 

321 

363 

34 

53 

73 

17 

261 

259 

311 

326 

364 

35 

54 

74 

18 

262 

260 

312 

365 

36 

55 

19 

263 

267 

313 

56 

20 

264 

268 

314 

57 

21 

269 

327 

(Cont . ) 


S85-1  1984  DENSITY  DATA  BASE  TAPES 
NINE  TRACK  6250  BP I  UNLABELLED  (W-I) 
(USE  SKIPF  OR  COPYBF) 


TAPE/FILE 

DAYS 

TAPE/FILE 

CC0290/1 

210,211 

CC0887/1 

/2 

212,213 

/  2 

/3 

214,215 

/  3 

/4 

213 

/  4 

/  5 

214 

/  5 

/  6 

215,216 

/  6 

n 

213 

/7 

/8 

214 

/8 

/  9 

216, 217 

/  9 

/10 

218, 219 

/10 

/II 

CC0828/1 

229,230 

/12 

/2 

231 

/ 13 

/3 

231,232 

/ 14 

/  4 

233 

/15 

/  5 

233,234 

/ 1 6 

/  6 

235 

/ 17 

/7 

235,236 

/18 

/  8 

237 

/9 

237,238 

/10 

239 

/II 

201 

DAYS 

TAPE /FILE 

DAYS 

239,240 

CC0517/1 

259,  260 

241,242 

/  2 

261 

241,242 

/  3 

261, 262 

243,244 

/  4 

263, 264 

245 

/  5 

265 

245,246 

/  6 

265,266 

247,248 

/  7 

267, 268 

246, 247 

/  8 

269,  270 

250, 251 

/  9 

219, 220 

252 

/10 

221 

247,248 

/II 

221,222 

249,250 

/ 12 

222,223 

251,252 

/ 13 

223,224 

253,254 

/ 14 

225,226 

254,255 

/15 

225, 226 

256,257 

/ 1 6 

227 

256,257 

/ 17 

226,227 

258,259 

/ 1 8 

228,229 

CC0833/1 

269,270 

/  2 

271,272 

/3 

272,273 

/  4 

274,275 

/  5 

276 

/  6 

275,276 

/  7 

277 

/  8 

280,281 

/  9 

282 

/10 

277,278 

/II 

279, 280 

/ 1 2 

271 

/ 13 

275 

/ 1 4 

281,282 

/ 1 5 

283 

y,\ 


ACCELEROMETER 


RAW/FILTER  TAPE  FORMAT 


HEADER  DESCRIPTION  F.QRHAT 

0.1  word  count  (40)  I 

0.2  group  count  (1)  I 

1  satellite  ID  A 

2  year  of  data  (since  1900)  (YY)  F 

3  Julian  date  of  data  (YYDDD)  F 

4  vacant  F 

5  order  of  temperature  fit  polynomial  F 

6  A0  (temperature  fit  polynomial  coefficients  )  F 

7  A1  F 

8  A2  F 

9  A3  F 

10  A4  F 

11  number  of  temperature  points  used  in  fit  F 

12  start  time  of  temperature  fit  F 

13  stop  time  of  temperature  fit  F 

14  number  of  missing  data  frames  F 

15  frame  increment  (time  between  xyz  samples)  F 

16  date  of  test  DMA  (raw  data)  run  (MM/DD/YY)  A 

17  Julian  date  of  raw  data  run  (YYDDD)  R 

18  start  time  of  accelerometer  data  F 

19  stop  time  of  accelerometer  data  F 

20  vacant  F 

21  start  time  of  filtered  data  F 

22  stop  time  of  filtered  data  F 

23  Julian  date  of  filter  run  (YYDDD)  R 

24  filter  length  F 

25  F3  X-axis  filter  parameters  F 

2  6  F4  F 

27  F3  Y-axis  filter  parameters  F 

28  F  4  "  F 

29  F3  Z-axis  filter  parameters  F 

30  F  4  "  F 

31  flag  for  treatment  of  missing  points  F 

32  number  of  missing  points  F 

33  flag  for  treatment  of  wild  points  F 

34  number  of  X-axis  wild  points  F 

35  number  of  Y-axis  wild  points  F 

36  number  of  Z-axis  wild  points  F 

37-40  vacant  F 
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ACCELEROMETER 


RAW/FILTER  TAPE  FORMAT 


DATA 

0.1 

0.2 

1 

2 

3 

4 

5 

6 

7 

8 

9-1024 


DESCRIPTION  FORMAT 

word  count  (8)  I 
group  count  (128)  I 
Greenwich  time  (seconds)  F 
X-axis  acceleration  (raw)  g's  F 
Y-axis  acceleration  (raw)  g's  F 
Z-axis  acceleration  (raw)  g's  F 
X-axis  acceleration  (filtered)  g's  F 
Y-axis  acceleration  (filtered)  g's  F 
Z-axis  acceleration  (filtered)  g's  F 
temperature  F 
cyclic  repetition  of  1-8  F 


ACCELEROMETER 


MERGE /EPHEMER IS  TAPE  FORMAT 


HEALER 

Same  as  raw/filter  tape 


'  ...■DATA  DESCRIPTION  FORMAT 

0.1  word  count  (30)  I 

0.2  group  count  (64)  I 

1  Julian  date  of  data  (YYDDD)  F 

2  calendar  month  F 

3  calendar  day  F 

4  Greenwich  time  (seconds)  F 

5  altitude  (Km)  F 

6  latitude  (  +  90  )  F 

7  longitude  (+E)  F 

8  geocentric  velocity  (Km/sec)  F 

9  velocity  relative  to  atmosphere  (Km/sec)  F 

,  10  local  time  (seconds)  F 

11  revolution  number  F 

12  angle  between  geocentric  and  atmospheric 

velocities  (degrees)  F 

13  vacant  F 

14  vacant  F 

,  15  solar  declination  (degrees)  F 

*  lr  vacant  F 

I  17  vacant  F 

18  vacant  F 

19  X-axis  acceleration  (raw)  F 

20  Y-axis  acceleration  (raw)  F 

21  Z-axis  acceleration  (raw)  F 

I  22  X-axis  acceleration  (filtered)  F 

23  Y-axis  acceleration  (filtered)  F 

24  Z-axis  acceleration  (filtered)  F 

25  temperature  F 

26  vacant  F 

27  vacant  F 

28  X-direction  clockwise  from  north  (degrees)  F 

29  magnetic  local  time  (seconds)  F 

30  vacant  F 

31-1920  cyclic  repetition  of  1-30  F 
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ACCELEROMETER 


MODEL/DENSITY  TAPE  FORMAT 


HEADER 

Same  format  as  raw/filter  tape,  except  model  run  date 
(MM/DD/YY)  added  in  word  37  and  ref.  area  M  added  in 
word  40 


DATA 

DESCRIPTION 

FORMAT 

0.1 

word  count  (30) 

I 

0.2 

group  count  (64) 

I 

1 

X-axis  wind  speed  (Km/sec) 

F 

2 

Y-axis  wind  speed  (Km/sec) 

F 

3 

Z-axis  wind  speed  (Km/sec) 

F 

4 

Greenwich  time  (seconds) 

F 

5 

altitude  (Km) 

F 

6 

latitude  (  +  90  ) 

F 

7 

longitude  (+e) 

F 

8 

geocentric  velocity  (Km/sec) 

F 

9 

atmospheric  velocity  (Km/sec) 

F 

10 

local  time  (sec) 

F 

11 

revolution  number 

F 

12 

angle  (in  degrees)  between  (8)  and  (9) 

F 

13 

CDX 

F 

14 

CDY 

F 

15 

normalized  density  (200  Km) 

F 

16 

J7 1  model  density 

F 

17 

MS IS  model  density 

F 

18 

measured  density 

F 

19 

X-axis  acceleration  (raw) 

F 

20 

Y-axis  acceleration  (raw) 

F 

21 

Z-axis  acceleration  (raw) 

F 

22 

X-axis  acceleration  (filtered) 

F 

23 

Y-axis  acceleration  (filtered) 

F 

24 

Z-axis  acceleration  (filtered) 

F 

25 

temperature 

F 

26 

Kp  (6-hour  lag) 

F 

27 

F  (one  day  lag) 

F 

28 

X-direction  clockwise  from  north 

F 

29 

magnetic  local  time  (seconds) 

F 

30 

CDZ 

F 

-1920 

cyclic  repetition  of  1-30 

F 
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